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Q\ ' Abstract 

cN : 

, We study the behavior near the singularity t=0 of Gowdy met- 

l/"") ■ rics. We prove existence of an open dense set of boundary points near 

which the solution is smoothly "asymptotically velocity term domi- 
nated" (AVTD). We show that the set of AVTD solutions satisfying a 
uniformity condition is open in the set of all solutions. We analyse in 
detail the asymptotic behavior of "power law" solutions at the (hith- 
i | erto unchartered) points at which the asymptotic velocity equals zero 

Jh . or one. Several other related results are established. 
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1 Introduction 

The Gowdy family of space-times [8] constitutes an interesting toy model to 
study formation of singularities in general relativity. This family of metrics 
is sufficiently simple to hope to analyse the resulting singularities in an 
exhaustive way. It is sufficiently non-trivial so that the relevant dynamical 
behavior has not been understood so far. The main question of interest 
is the curvature blow-up - or lack thereof - at the boundary t = of the 
associated space-time. The reader is referred to [6] for a further discussion 
of this issue, we simply note that the relevant geometric information can be 
obtained by deriving a sharp asymptotic expansion of the solutions near the 
singular set t = 0. The main purpose of this work is to prove a stability 
result for the existence of such expansions. 

In Gowdy space times the essential part of the Einstein equations reduces 
to a nonlinear wave-map-type system of equations [8] for a map x from 
(M,g a/3 ) to the hyperbolic plane (J$?,h ab ), where M = [T, 0) x S 1 with 
the flat metric g = —dt 2 + d9 2 . The solutions are critical points of the 
Lagrangean 

^\A = \l tg a ^h ab d a x a dpx b d9dt . (1.1) 
2 J M 

This differs from the usual wave-map Lagragean by a supplementary mul- 
tiplicative factor t. It is sometimes convenient to use coordinates P,QgR 
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on the hyperbolic plane in which the hyperbolic metric h a b takes the form 

h = dP 2 + e 2P dQ 2 . (1.2) 

Let X t = Qft,Xg = D denote the Levi-Civita connection of h a b , and 
Dg = := Dx g , D t = jji := Dx t - The Euler-Lagrange equations for (1.1) 
take the form 

DX t DX e _ X t 

~T>T ~d¥ ~ T [L6) 

or, in coordinates, 

ban r. 9fX a 



a x a + Tl c o xd^d^x 



t 



where the T's are the Christoffel symbols of h a b, and □ = df — di. Global 
existence of smooth solutions on (— oo,0) of the Cauchy problem for (1.3) 
has been established by V. Moncrief [13]. 

For further use we note the non-vanishing Christoffel symbols of h: 

^ Q = -e 2P , rgg = rg P = i. (1.4) 

In the (P, Q) coordinates one thus has 

d 2 P -d 2 P=-^f + e 2P ((d t Q) 2 - (dgQf) , 
d 2 Q -d 2 Q = -^--2 (d t Pd t Q - dgPdgQ) . 

We consider solutions defined on sets Cl(a, b, to), where 

t < , a < b , U(a, b, t ) := {t < t < , a + t < 9 < b - t} (1.5) 

(see Figure 1). Thus our analysis is local, if the solution satisfies certain 
properties on an interval [a, b] C S 1 , then the conclusions hold on that 
interval. Throughout this work we assume that the initial data for the map 
x at to are smooth functions of 9. We prove the following theorem (the 
Geroch group is defined in Section 4; the position function ip^ is defined by 
Equation (3.24)): 

Theorem 1.1 Let (x(to, ■), Xt(to, ■)) be Cauchy data for a solution of the 
Gowdy equations on £l(a,b,to) such that the associated solution x has uni- 
formly controlled blow-up; by this we mean that 



sup 

9e[a-\t\,b+\t\] 



(j2 \t k+1 D k gXg\ + \t k+l D k e X t \\ (t, 9) , (1.6) 

\fc=0 k=l / 



with 1 = 2. Then: 
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Figure 1: The set Q(a,b,to). 



(i) For all 9 E [a, b] the function \tXt\(t, 9) converges to a velocity function 
v(0) as t tends to zero, uniformly in 9. 

(ii) There exists an open dense set on which v is smooth. 

(iii) [a, b] can be covered by a finite number of intervals [aj,6j] with the 
following property: for each i there exists an element Gi of the Ge- 
roch group such that GiX has a smooth velocity function < v < 1 
and a smooth position function ip^, except perhaps on the boundary 
of the set {v(9) = 0}. Further GiX satisfies a power law blow-up, 
Equation (8.20). 

(iv) In the associated space-time the curvature scalar R a p^^R a ^ & blows up 
in finite proper time on every causal curve approaching 

08 := {0} x ([a, b] \ {v(9) = 1}) x^xS 1 . 

In particular the associated Gowdy space-time is inextendible across 

m. 

(v) There exists r] > such that for all initial data (x(tQ, -),X t (to, •)) sa ^' 
isfying 1 

\\(x(t ,-) - x(t , -),X t (t , •) - X t (t ,-))\\ H 3 m 2 < T) 

the associated solution x also satisfies (1.6), and hence also (i)-(iv) 
above. 

x We equip fi(a, 6, to) with the Riemannian metric dt 2 + dd 2 , this together with the 
metric h on Jffi induces Riemannian metrics on all the bundles involved. We use those 
metrics to measure the distance between points on those bundles. 
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Remark 1.2 Actually it suffices to have a sequence U — ► along which (1.6) 
holds. A recent result of Ringstrom [18] 2 can be used to lower to £ = 1 the 
threshold I in (1.6). 

The proof of Theorem 1.1 can be found at the end of Section 12. 

It is of interest to enquire whether the uniform blow-up condition (1.6) 
is necessary for AVTD behavior of the solutions. Consider, for example, 
an AVTDg^'^ solution, as defined in Section 3, for which the error terms 
in (3.5)-(3.6) and in their derivative counterparts (see (3.10)) are uniform 
in 8. If v\ is strictly smaller than one (no negative lower bound assumed), 
then the solution satisfies (1.6) 3 . This shows in particular that the set of 
solutions satisfying the hypotheses of Theorem 1.1 is not empty, as existence 
of a large class of AVTD^'^ solutions satisfying v\ < 1 follows from the 
results in [14]. 

The second main result of this work is the proof that for every solution 
there exists an open dense set Cl C S 1 near which we have complete control 
of the solution: 

Theorem 1.3 Consider a solution x defined on Q(a,b,to). There exists an 
open dense set Q C [a, b] such that x is in a neighborhood of 

{o} x n. 

The examples discussed in Section 3 show that the result is sharp, with 
the following proviso: the known examples have a velocity function defined 
everywhere, even at points where it is not continuous, while Theorems 1.3 
and 12.1 leave open the possibility of existence of points where the velocity 
is not defined. Such points are characterised in point (i) of Proposition 12.3. 

The third main result of this paper is an exhaustive analysis of the 
asymptotic behavior of power-law solutions at points at which v vanishes, 
or equals one. This last case is especially important for the discussion of 
strong cosmic censorship, we refer the reader to [6] for applications. We 
note that no results concerning those velocities were available so far in the 
non-polarised case. 

2 The results presented in [18] have been made available as a preprint [17] while this 
paper was being prepared for publication. 

3 This follows immediately from Proposition 1.4 together with the calculations of the 
proof of Lemma 8.7. For points 6 at which v\ > 1 one expects Qoo to have vanishing 
derivatives, compare point (ii) of Proposition 3.1 and Theorem 10.2. Then, if vi{9o) > 1 
and if <9gQoo(#o) = 0, i = 1,2,3, then (3.5)-(3.6) and their derivative counterparts give 
pointwise decay of the function appearing under the sup at the left-hand-side of (1.6), but 
uniformity is far from being clear. 
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The results discussed above are established through a series of auxiliary- 
results which have some interest in their own. We say that a solution x sat- 
isfies a power law blow-up, or is of power-law type, or is a power-law solution, 
if the norm of the theta derivatives vector \Xg\ does not blow up faster than 
for some positive constant e, when approaching the singularity t = 0. 
All solutions of the smooth Cauchy problem on T 3 analysed in detail so far 
satisfy 4 a power law decay. 

It is simple to show that every solution with a power law decay has a 
continuous asymptotic velocity function v (see the proof of Theorem 8.3 
below). The associated solutions of the vacuum Einstein equations have 
curvature blowing up uniformly, except perhaps at the set of points 9 at 
which v{9) = 1. Consider the set of initial data for solutions satisfying 
a power law decay and for which v < 1, uniformly in 9. We show - see 
Theorem 11.1 below - that this set is open in the set of all initial data; this 
is one of the steps of the proof of Theorem 1.1. We further show that for 
those solutions v is smooth except perhaps at the boundary of the set of 
points at which v vanishes. Theorem 11.1 leads to a sharper version of the 
stability of the singularity theorem for (2/3,2/3,-1/3) Kasner metrics, see 
Theorem 9.1. 

An important element of our analysis is the action of the Geroch group, 
as defined in Section 4. In fact, the key ingredients of our analysis are the 
results in [4] together with the following: 

(i) The analysis of the action of the Geroch group in the work of Rendall 
and Weaver [15]; 

(ii) The reformulation of the wave-map equations as a first-order system 
of scalar equations by Christodoulou and Tahvildar-Zadeh [3]; 

(iii) The small-derivatives stability result of Ringstrom [16]. 

We finish this introduction by recalling some results from [4] which will 
be useful in the sequel: 

Proposition 1.4 (Time-weighted pointwise estimates; Proposition 3.2.1 in [4]) 
Let x(t ,9) G C fe (5 1 ), k > I, X t (t ,9) e C fc " 1 (5 1 ). For all t > t we have 

4 More precisely, the examples known to us satisfy a power law decay on all sets Ct (ip) 
as defined in Equation (2.3) below. The constants are uniform in tfi away from the points 
at which the asymptotic velocity has spikes, or discontinuities, or crosses zero or one. 
In [6] an explicit self-similar solution has been given which does not satisfy the power law 
decay, but this solution does not fit into a Cauchy problem framework. 
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(i) (|^ t | 2 +|^| 2 )(t,^) < 2{sup^ e[£ ,_, + , 0> ^_, o] (|^ t | 2 + |^| 2 )(i ,^)} (^) 2 . 

(ii) If k > 2, then there exist constants C depending only upon the argu- 
ments listed such that, for all 1 < \a\ < k, 

|^ a x|(t,e)<C(|a|,*oJ|X e (to)|| cN -iJ|X t (*o)|| c |a|-i)|tr |Q| . (1.7) 

Remark 1.5 It has been pointed out to us by H. Ringstrom that the proof 
of Proposition 3.2.1 in [4] (compare [4, Equation (3.2.5)] together with the 
argument leading to Equation (3.2.9) there) actually gives an inequality 
somewhat stronger than (i) above: 

£j 2 (\x t \ 2 + \x e \ 2 )(t,e)< 

-{ sup \x t -X e \ 2 (t ,ij)+ sup \X t + X e \ 2 (t ,i;)} . 
1 y 4>e[e-t+t ,e+t-t } ^e[e-t+t ,e+t-t ] J 

(1.8) 

Equation (1.8) carries more information about the solution than (1.7), which 
can be seen e.g. when the initial data have small ^-derivatives. 

PROPOSITION 1.6 (Time- weighted Sobolev decay; Proposition 3.3.1 in [4]) Let 
x e C i ([t ,0) x S 1 ) and let X e (t ,-), X t (t ,-) G H^S 1 ), i > 1. Then there 
exist constants depending only upon the arguments listed such that 

(i) For all 1 < |a| < i + 1, 

g^(t) = j> de\t\ 2 W\D a x\ 2 < C(\a\, ||X*(to)||ff |a| _ l( si), ||X t (to)||ff |a| _ l( si), t ). 

(ii) // at least one differentiation is a 6 differentiation we have 

lim g (a) (t) = 0. 

(hi) If at least one differentiation is a 9 differentiation then 9 ^ G i 1 ([to> 0]) 
and 

f ^rfi ds ^ C'(\a\,to, \\Xe(t )\\ HH _ l{s ^ ||^(to)|| ff|a| _ l( si))- 

Jto \ s \ 
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Figure 2: The truncated domains of dependence C\ {tp). 

2 Problems with 6 derivatives, self-similar solu- 
tions 

As already mentioned in the introduction, all published solutions of the 
Gowdy equations known to us, and for which the asymptotic behavior is 
reasonably well understood [4, 5, 10, 14-16], have the property that 

\tXg\ < Ct e (2.1) 

for some e > 0, with the bound being optimal 4 . The power law is very 
useful for the control of the analytic properties of the solutions, but it is 
not necessary for curvature blow-up. In any case the bound (2.1) certainly 
implies, for all ip G S 1 , 

[ \X e \ 2 dtdO <oo , (2.2) 

where for to < the set C t ° (tp) is defined as (compare Figure 2) 

C t ° (V) = {to < t < , < ~ V> < \t\} ■ (2.3) 
We shall say that lim^o = a if 

lim sup \f(t,9)-a\=0. (2.4) 

-\t\<0-tf><\t\ 

Such limits look a little awkward at first sight; however, they arise naturally 
when considering the behavior of the geometry along causal curves with 
endpoints on the boundary t = 0. Further, such limits appear naturally in 
our results below. 

We have a partial converse to (2.1): 
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PROPOSITION 2.1 (Proposition 3.4.1 in [4]) At every ip at which (2.2) holds 
we have, for all multi-indices a, 

lim C oW \t\^ +1 \D a X d \=0. 

Remark 2.2 In Section 6 below we give further integral conditions which 
ensure pointwise convergence of \tX t \ to a number v{ip). Yet another crite- 
rion for existence of v(ifj) is given by Proposition 12.3. 

Proposition 2.1 begs the question of existence of solutions for which (2.2) 
fails. An obvious candidate is given by self-similar solutions: 

x(t,0) = y(0/t), (2.5) 

for some map y from M to the hyperbolic space. It would be of interest to 
find all solutions satisfying (2.5). Here we note the following family of such 
solutions: let a, f3 G M and let V : R — > J$? be an affinely parameterised 
unit-speed geodesic in Jif, for \9\ < —t set 

x(t, 9) = r(a arcsin (j) + /?) • ( 2 -6) 

It is easily checked that (2.6) solves the Gowdy equation (1.3). Equa- 
tion (2.2) fails for the solution (2.6) when a / 0, as expected. It turns 
out that the solutions (2.6) do not fit into the Cauchy problem framework 
because they are singular on the whole light cone \9\ = —t, while the so- 
lutions we are interested in are smooth at t = to- Singular solutions can 
sometimes be used to produce smooth examples of bad behavior, but we 
have not managed to exploit this solution to do that. In view of our stabil- 
ity results here it would be important to construct a solution of the Cauchy 
problem where (2.2) fails, or to prove that such solutions do not exist. 

3 AVTDfc behavior, spikes and discontinuities 

In [5, 11, 14] a large class of solutions of (1.3) has been constructed with the 
following behavior: 

P(t,9) = -vi(0)ln|*l + ^°o(0)+o(l) , 0<v 1 (0)<l, (3.1) 
Q(t,e) = Qoo(e) + \t\^W(l, Q (9)+o(lj) . (3.2) 
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A generalisation of those formulae to arbitrary velocities requires a careful 
study of the field equations. For instance, an analysis of the indicial ex- 
ponents of the linearised equations suggests the following behavior of Q at 
points where v(9) = n £ N* (compare [6] for n = l) 5 



f2n 

Q(t,6) = Q^(e) + —dl n Qoo(e)ln\t\+iP Q (6)t 2n + o(t 2n ). (3.3) 

While there is no existence statement for solutions with a non-zero coeffi- 
cient in the In \t\ term above, we expect that such solutions can actually be 
constructed. We note that if v G N* on an interval, then no log term will 
occur on that interval. 

Applying the "solution-generating transformation" (3.19) described be- 
low to (3.3) at a point at which one further has Qoo(0) = leads to a solution 
(P',Q') with a negative P- velocity v[(9) = —n and with a logarithmically 
blowing-up Q function 

Q'(t,e) = e 2P ^{j^dl n Q^{e)\n\t\+^ Q (e)^+o{l). (3.4) 

The above discussion suggests that the following will capture the asymptotic 
behavior of a large class of solutions of the Gowdy equations: 

P(t,9) = -v 1 {6)\n.\t\ + P 00 (9) + o(l), (3.5) 

\t\ 2 ^U Q (9) + o{\)) , o<^)p ; 

|tr iW (<2ln(0) In \t\ + 4> Q {e) + o(l)) , < v\(9) G N ; 

Qin(0)ln|t| + o(l), «i(0)€-N*; 
o(l) , -N* ^ vi(9) < . 

(3.6) 



Q(t,9) = Qoo(9) + < 



The function 



v := (3-7) 



will be called the velocity function, while Qoo will be called the Q -position 
function. Those functions have the following geometric interpretation [9]: 
for v > the path 

r - T (t) := (P(- e - T ,9),Q(-e- T ,9)) 



5 The expected a priori estimate \t 1+1 D^Xo] — > implies that d' g Qoo(0) — for i 
1, . . . , 2n — 1, which is implicit in (3.3). 
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approaches - in a sense made precise by (3.5)-(3.6) - the affinely parame- 
terised /i-geodesic 



T^T e {r) -(viWr-tpiO^QooiO)) , 

with v - the length of the velocity vector of Tg. The point Qoo(6) is then 
the uniquely defined point on the conformal boundary of the hyperbolic 
space at which the geodesic Tg accumulates. Clearly this interpretation 
breaks down at v(9) = 0, which suggests that solutions might display strange 
features, not necessarily compatible with (3.5)-(3.6), at the boundary of the 
set {v{6) = 0}. 

We shall say that a map x = (P, Q) is in the AVTD^'^) class if there 
exist functions v±, Qoo, and Q\ n such that 

P(f,9) = -«i(0)ln|t| + , ( 3 - 8 ) 

C(t9) f Qoo(9) + o(l) , v^-W; 

W ' J 1 Q\n(Q) In \t\ + Qoo(9) + o(l) , vi <E -N* . [6 - y> 

We shall say that a solution is in the AVTD^ P ' < ^' ) class if (3.5)-(3.6) hold with 
functions v\ , P^ , Qoo , Q\ n and i/jq which are of differentiability class (on 
closed intervals the derivatives are understood as one-sided ones at the end 
points). For the purposes of the AVTDj.^'^ definition the function Q\ n is 
assumed to be extended by to the set v±(9) — N*; we emphasise that 
such an extension will not be assumed in Definition 3.4 below. For k > we 
will assume that the behavior (3.5)-(3.6) is preserved under differentiation 
in the following way: 

V0<i+j<k &> 9 {td t ) l (p{t,9) + v l {9)\n\t\ -Poo(0)) =o(l) , (3.10) 

similarly for Q. 

Note that the classes AVTD( P ' Q ) and AvTd{, P ' Q) do not coincide. 

Unless explicitly stated otherwise the o(l) symbol denotes pointwise con- 
vergence to zero as t tends to zero at fixed 9. Similarly O(l) means uniformly 
bounded in t at fixed 9. An alternative meaning of o(t), which will be some- 
times used, is provided by uniform convergence to zero on the set 0(a, b, to) 
as defined by (1.5). Yet another possibility is convergence to zero in a 
lim c o (y,) sense, as defined by (2.4); in any case we will make precise state- 
ments when needed. 

Using the above solutions, Rendall and Weaver [15] have constructed so- 
lutions of (1.3) which display "spiky features" . They separate those solutions 
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into two classes, one called "false spikes" and one called "true spikes". An 
instructive example of this behavior is provided by a family of Gowdy maps 
discovered by Moncrief [12], and analysed in detail in [4, Appendix B]. They 
are given there in terms of the polar coordinates on the hyperbolic space, 

h = dp 2 + sinh 2 p dip 2 , (3.11) 

which are related to the (P, Q) coordinates by the formulae 6 

e p = cosh p + sinh p cos <p 

= ^(e^l+cos^+e^l-cos^)) , (3.12) 
e p Q = sinhp simp . (3.13) 

Moncrief's ansatz 

p = p(t), tp = n6, n G N , (3.14) 

leads to the following: every solution is uniquely determined by two numbers 
e [0, 1), poo € R, such that 

p = - Voo ln\t\ + Poo + o(l) . (3.15) 

Inserting (3.14)-(3.15) into (3.12) one finds 



P(t,9) 



{ -«oolii|t| +Poo+ln( 1+c ° 2 s( " e) ) +o(l) , uO^tt mod27r; 
[ Voo hi \t\ — poo + o(l) , n9 = 7T mod 2tt, 

(3.16) 

Qm = JrSS)+°( 1 )' mod2vr; 

1 , n9 = ir mod 2-7T. 

If we define «i : S 1 — > R by the equation 

«i(0) :=lim|t|P t , (3.18) 

then = > except at n isolated spike points 9 m = (1 + 2m)7r/n, 
m € Nn [0,n], at which fi is equal to — Voo < 0. Equation (3.16) shows that 
the subleading term in P blows up logarithmically at the spike points, so 
that no uniformity in 9 for that term can be expected near those points in 
norms which control pointwise behavior of P and Q. It is interesting that 

We are very grateful to Marsha Weaver for several enlightening discussions concerning 
the issues discussed in this section, and for providing us formulae (3.12)-(3.13). 
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even though blows up as one approaches the spike points, it is finite, 
actually vanishing, there. Next, even though at each fixed 9 we have 

lim P e (t, 9) = , 

there are timelike curves reaching the boundary along which Pq does not go 
to zero: for example, for n = 1, 

\imP e (t,7r + ae- p M) = oo , lim P e (t, ir + ae" 2 ^) = | . 

The above examples provide solutions with an arbitrary finite number of 
spikes. Solutions with a countably infinite number of spikes accumulating at 
some point 9^ £ S 1 can be constructed as follows: Consider an AVTD^ > ' ( ^" ) 
solution such that the function Qoo in (3-6) has an infinite number of isolated 
zeros Oi accumulating at 9^, and such that v\ avoids zero in a neighborhood 
of #oo; the existence of such solutions follows from [14]. Following [15], one 
then performs the following "inversion" of the hyperbolic plane: 

e P = Q2 + e -2P > Q' = Q2 + £ -2P ■ ( 3 - 19 ) 

This is an isometry of h and therefore maps solutions into solutions. It is 
easily seen that (P',Q r ) will have a spike at each of the points 9i, yielding 
the desired spiky solution. 

The discontinuities discussed so far consisted of isolated points at which 
vi changes sign. Solutions with jumps of v \ can be constructed as follows: 
consider any AVTB^'^ solution such that the zero set of the function Qoo 
in (3.6) is a closed interval [a, b], with v\ strictly positive near the end 
points. As before, the existence of such solutions follows from [14]. It is 
easily seen from (3.19) that the velocity function v[ associated with the map 
x' = (P',Q') will jump from v\(a) to — v\(a) at a, and will be continuous 
from the right there. Further, as 9 increases from a to b the new velocity 
function v[ will continuously attain the value — v\{b) when b is approached 
from the left, and jump to v\{b) immediately afterwards. 

Clearly, the above behaviors can be combined to give infinite sequences 
of pointwise jumps and/or intervals on which v[ is negative, with the set of 
discontinuities of v[ accumulating at a given point. 

In fact, let F C S 1 or F C [a, b] be any non-empty closed set without 
interior, we claim that there exists a smooth function ipp such that 

^({O}) = F . 
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In order to see this, let x £ F, and let (x-,x+) be the largest open interval 
containing x which does not meet F (hence x± € F), we set a(x) = (x — 
x_)(x + — x). Define 

<Pf(x) = { e li M x) ) otherwise . 

Then 93^ has all the required properties. Using the function (fp as Qoo, with 
ui equal, e.g., to the constant function 1/2, after performing an inversion 
we obtain a new function v[ which equals —1/2 on F, and 1/2 on S 1 \ F. If 
F is a fat Cantor set one obtains a rather wild set of spikes, with measure 
as close as desired to that of S 1 or that of [a, b] by choosing F suitably. 

The spikes discussed so far are called false spikes, as they can be thought 
of as an artifact of the (P, Q) coordinate system above: no discontinuous 
behavior occurs in the (p, tp) representation of the solutions. 7 However, the 
(P, Q) coordinates are very useful when analysing the Gowdy-to-Ernst trans- 
formation, because that transformation has a very simple form precisely in 
the (P, Q) coordinates: given a solution x = (P, Q) of the Gowdy equa- 
tions, one defines a new solution x by performing the "Gowdy-to-Ernst" 
transformation [15]: 

P--P- In |t| , e p 8 t Q := -e p d e Q , e p d e Q := -e p 8 t Q . (3.20) 

The new map satisfies again the Gowdy equation (1.3). As shown by Rendall 
and Weaver, this has significant consequences: By definition, a true spike is 
the image of a false spike after a Gowdy-to-Ernst transformation has been 
performed; equation (3.20) shows that any discontinuity in v\ leads to a 
discontinuity in the velocity v\ associated with the map (P,Q). For AVTD 
solutions this typically leads to a discontinuity in the geometric velocity 
function v = \v\\. For instance, for the solutions (3.14) the transformation 

7 Applying isometries of the hyperbolic plane to the solution has the effect of reshuffling 
Killing vectors, and can thus be considered as an irrelevant "coordinate transformation" 
from the space-time point of view. We note that the isometry (3.19) changes the orienta- 
tion of the hyperbolic plane. The accompanying relabeling of the Killing vectors changes 
the space-time orientation, leading thus to a non-equivalent solution if a space-time ori- 
entation has been chosen. However, we can always perform a second inversion about a 
different point, regaining the original orientation. If the map 9 Q x (0) is surjective, 
this will always introduce at least one "false spike" in the transformed solution. In other 
words, in the surjective case there will be no (P, Q) representation of the solution without 
"false spikes" . This gives some geometric meaning to those. 
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(3.20) leads to 



p = \ -(1 - Voo) In 1*1 " Poo - In ( 1+co 2 s(wg) ) + o(l) , / vr mod 2vr; 
[ — (1 + Voq) In |f| + poo + o(l) , n6 = it mod 27T, 

(3.21) 

which clearly results in a \v\ which is not continuous at n6 = it mod 2ir. 

We refer the reader to [15, Section 6] for a further discussion of iterations 
of the above. 

It is interesting to enquire about independence of the conditions (3.8)- 
(3.9). It turns out that (3.8) is the key requirement, up to a Gowdy-to-Ernst 
transformation: 

PROPOSITION 3.1 (i) At each point 9 at which (3.8) holds with v\ > we 
also have (3.9). If the term 0(1) in (3.8) is uniform in 9 over some 
interval I then Qoo is continuous on I. If further v\ is uniformly 
bounded away from zero, then the term o(l) in (3.8) is uniform in 
9 el. 

(ii) If vi > 1 on an interval I, with the term 0(1) in (3.8) uniform in 9, 
then there exists a constant A such that 



V# €/ Q(t,9)^ A 
(we say that x is asymptotically polarised on I). 

(hi) At points 9 or intervals I on which (3.8) holds with v\ < 1 the conclu- 
sions of point (i) above hold for the Gowdy-to-Ernst transformed map 
x. If v\ < on I then the conclusion of point (ii) holds for Q. (In the 
case of an interval I we assume that the term O(l) in (3.8) is uniform 
in 9.) 

Proof: Point (i) of Proposition 1.4 shows that 

\Qt\ + \Qo\ < C{9 1 ( : P < C(9)\tr^' 1 . (3.22) 

Integrating in t one obtains (3.9). If C{9) can be made ^-independent, then 
Qoo is a uniform limit of continuous functions, and therefore continuous, 
which gives (i). If v\(Q) > 1, then (3.22) shows that \Qg\ tends to zero as t 
goes to zero, which easily implies that Qoo is constant over /. This proves 
(ii). Applying the Gowdy-to-Ernst transformation (3.20) one finds that the 
hatted velocity function v\ associated with P equals 1 — v\, and point (hi) 
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follows. 



□ 



Let us summarize the properties of the (P, Q) coordinates which follow 
from the above: 

(i) They provide a simple explicit formula for the Gowdy-to-Ernst trans- 
formation; 

(ii) They describe faithfully the geometric behavior of the solutions except 
near those points at the conformal boundary with <p = n mod 2tt; 

(hi) They describe faithfully the geometric behavior of those solutions on 
intervals of 9 on which v vanishes. 

The problem with ip = it mod 2ir above is avoided by turning to the (p, ip)- 
description of the solutions. So, instead of (3.5)-(3.6) we write 



(We impose the usual restriction that p > so v above is necessarily non- 
negative, though in some situations it might be convenient not to do this, 
allowing p to be negative, but then identifying the points (p, <p) with (— p, p— 
7r).) It follows from (3.12) that the notation v for the p- velocity in (3.23) 
is comptatible with (3.5) and (3.7). A map x = (p, p) will be said to be 
avtd(^) on an interval [a, b] if there exist real valued functions v and <poo 
such that for 9 G [a, b] we have 



x will be said to be AVTDjf' w on an interval [a, b] if (3.23)-(3.24) holds 
with functions v, p^, <poo, ip\ D and ip^ which are C k on [a, b}. For k > the 
derivatives are assumed to behave as in (3.10). 

Because of the usual polar coordinate singularity at p = the (p, p) 
coordinates do not always reflect the geometric character of the dynamics 
for solutions on intervals on which v(9) = 0. 



p(t,9) 




p(t,9) 



p(t,9) 
<p(t,0) 



-v(6)ln\t\ + 0(1) , 

Poo(9) + O(l) . 



(3.25) 
(3.26) 
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We note the following result, which follows immediately from the calcula- 
tions in the proof of Proposition 4.2 below. It shows that the only points 9 at 
which the distinction between AVTD^'^ and AVTDj^'^ behavior matters 
are those at which v vanishes (where (p, ip) might be singular) or at which 
(poo = tt mod 2tt (where the restriction of Q to the conformal boundary of 
the hyperbolic space is singular): 

Proposition 3.2 Let k > 0. 

(i) If the map x is AVTD^'^ on [a\, b\] with v\ avoiding zero on [a±, b\\, 

then it is AVTD^' V \ Further ip^ mod 2tt avoids tt on [ai, b\] if v\ > 
0, while p^ = tt mod 2tt if vi < 0. 

(ii) If the map x is AVTD^'^ on [ai,b\] with ipoo mod 2ir avoiding tt on 
[a x M\, then it is AVTD { k P,Q) . 

Remark 3.3 We have an obvious equivalent of points (i) and (ii) of Propo- 
sition 3.1 for the (p, <p) representation of the solutions, with identical proof, 
regardless of whether or not (p^ meets tt mod 2tt. 

What has been said so far in this section leads naturally to the following 
definition: 

Definition 3.4 A map x : ft(a, b, t ) — ► J#2 will be said to be AVTD, respec- 
tively AVTDk, on [ai,bi] C [a, b] if there exists a function v : [ai,i>i] — > M + 
such that: 

(i) For every interval I C [ai,6i] on which v vanishes the map x is 
AVTD^ Q \ respectively AVTD { k ' Q} , near {0} x /, with \v±\ = v. 

(ii) For every interval I C [ai,6i] on which v has no zeros the map x is 
AVTD^'f), respectively AVTD^ ,(p \ near {0} x /. 

4 The Geroch group and its action 

We consider solutions of the Gowdy equations defined on tt(a, b, to) (see 
(1.5)), for some a < b, to < 0. We fix 8 once and for all a (P,Q) coordinate 
system on (J#2,h). Following [7], we define the Geroch group W as the 
set of finite strings of the form G = G1G2 ■ ■ ■ Gfe, where each of the Gj's 

We emphasise that there are several coordinate systems (P, Q) in which the metric h 
takes the form (1.2), differing from each other by an isometry of h. 
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is either an isometry of (J#2,h), denoted by ij, or is a Gowdy-to-Ernst 
transformation (3.20), denoted by E. The Geroch group acts on solutions 
as follows: First, G\ - • • G n acts on x by first acting with Gn on x, then 
acting with G n ~\ on G n x, etc. Next, if x is a solution of (1.3), then we 
start by writing it in the (P, Q) coordinate system just chosen. Isometries 
act on solutions by composition. This implies that the ^-group product 
J1J2 of two isometries I\ and I2 is the composition I\ o I 2 of I\ with I2. 
The action of a Gowdy-to-Ernst transformation £ona solution is defined 
as follows: we integrate (3.20) with the integration constant chosen so that 
Q(to, a) = Q(to, a). This leads to the group product E 2 = ld^ 2 , where Id^ 2 
is the identity isometry of the hyperbolic plane. The £f -group products EI 
and IE, are defined by the above action on solutions. 

Let G = G1G2 ■ ■ ■ Gk, then any two adjacent Gowdy-to-Ernst trans- 
formations can be canceled out, leading to a shorter presentation of G. 
Similarly any two adjacent isometries can be replaced by a single isome- 
try. This leads eventually to a presentation of G such that isometries and 
Gowdy-to-Ernst transformations alternate. The number of Gowdy-to-Ernst 
transformations in the resulting presentation of G will be called the order 
of G. Thus, the order of / is zero, the order of E, or IE, or EI, is one, etc. 

The behavior near t = of all the known to us solutions of the Cauchy 
problem for the Gowdy equations is captured in the definition of the set 
below; the set ^2 is then a subset of ^1 with a genericity condition; we 
expect % to be useful in the analysis of the strong cosmic censorship prob- 
lem in the class of Gowdy space-times. The following comment is in order 
here: as emphasised in the previous section, the (P, Q) variables provide 
a parametrization of the hyperbolic space which does not always correctly 
reflect the geometric aspects of the asymptotic behavior of the solutions, so 
the reader might wonder why to invest so much effort to characterise the 
asymptotics of the Gowdy solutions in terms of those variables rather than, 
say, the {p, tp) variables of (3.11). The answer is that the Gowdy-to-Ernst 
transformation takes a simple form in the (P,Q) variables, compare (8.33), 
and this is what forces us to carry out the analysis below. 

Definition 4.1 1. Let %\ be the set of smooth solutions of the Gowdy equa- 
tions defined on Q(a,b,to) for which (3.5) -(3.6) holds with some functions 
v\, Poo, Qoo defined on [a,b], a function ipQ defined on the set {v± > 0}, 
and a function Q\ n defined on the set {v\ G satisfying the following: 

(i) v\ is uniformly bounded. 

(ii) v\ is continuous on an open dense subset of [a, b] . 
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(iii) The restrictions o/Qoo to the sets {9 : v±(9) — N} and {9 : v\(9) 6 
— N} are continuous functions on those sets. 

(iv) The restriction of P^ to the set {9 : Qoo and v\ are continuous at 9} 
is a continuous function on this set, similarly for Q\ n and ipQ. 

(v) At points at which Q\ n {9) / the error terms o(l) in (3.5) -(3.6) have 
the property that they remain o(l) after multiplication by \n\t\. 

2. We define % to be the subset of consisting of those solutions for 
which the sets of discontinuities and of critical points of v\ and of Qao are 
finite. 

It is an open question whether there exist solutions of the smooth Cauchy 
problem for the Gowdy equations which are not in 
In our analysis below we will need the following: 

Proposition 4.2 <ft\ and % are invariant under the action of isometries 
of(J%,h). 

Remark 4.3 The arguments of the proof of Theorem 12.1 show that the 
set of solutions which are in <2q is stable under the action of the whole 
Geroch group. It is not completely clear what happens with the action of 
the Geroch group on solutions in % , since the integration of the Q equations 
might introduce non-generic behavior. 

Proof: Let <j) be an isometry of the hyperbolic space into itself. We will 
write 

x = <p o x , 

and denote by (P, Q) the associated coordinate functions. 

If #o £ [a, b] is such that tq(6>o) = 0, then continuity of (p shows that (3.5)- 
(3.6) does hold for (P(t,9),Q(t,9)) at 9q, with the new velocity t>i(#o) = 
= vi{9). 

Let V be the open dense set on which v\ is continuous, then "V is a 
nonempty countable union of open intervals Ii,'f = Dili. By definition the 
velocity function v\ is continuous on each Ij. We rewrite the ij's as 

I t = t Vl > 0} U {vi < 0} U {v! = 0} . 

V v ' V v ' v ' 

h+ U- ho 

The simplest set to analyse is Ijo : f° r (t,9) E [to,0) x /jo the map x(t,9) 
stays in a compact set, so does <f> o x for any isometry 0, and the property 
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that (Poo, Qoo) are continuous on Jjo is clearly preserved under the action of 
isometries. 

Let us use the angle ip of (3.11) to parameterise the conformal boundary 
of the hyperbolic space, then to every point for which v\(9) > we can 
assign a unique (foo(9) such that the trajectory t — > x(t,6) asymptotes to 
the point at infinity ip oo (0). From (3.12)-(3.13) we have 

Vl{6) > =>- Qoo{0) = — -. -rrrr , (4.1) 

1 + COS(^oo(6>)) 

except at 

if ^(9) = ir mod 2vr (4.2) 

where a more careful analysis is required. We emphasise that a possible 
singularity arising here would only reflect the singular behavior of the Q- 
parameterisation of the conformal boundary, and not a singularity of Poo{9): 
continuity of ip^ on the set {v > 0} can be established as in Proposi- 
tion 3.1, working directly in the (p, ip) coordinates, regardless of whether or 
not cos(#) = — 1. In any case Qoo is continuous on ij + by hypothesis, and 
so is therefore <^oo- 

Inverting (3.12)-(3.13) one finds 



e" p + e p (l + Q 2 ) , (e- p + e p {l + Q*Y 2 



2 / + V 2 — 1 - 1 " ^ 



(The alternative solution 



always leads to a negative p. Equations (4.3)-(4.4) reflect the fact that a 
point (P, Q) corresponds both to (p, ip) and (— p, if — it). In the current proof 
we follow the usual convention that p > 0.) It follows that 

p(t, 6) = - Vl {9) In \t\ + Poo (9) + o(l) , (4.5) 

for some number p<x,(6). Inserting (4.5) into (3.13) one is then led to 



<p(t,0) = <p oo (0) + 



\t\ 2 ^ (%{9) + o{\)) , V!(9)^N*, 

\t\ 2vi{d) (<A„(0) In |*| + M°) + (!)) > «i W e ^* . 

(4.6) 
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for some numbers ip<p(0), </?in(#)- As tp^ is a continuous function of 9 on 
continuity in 9 of Poo(#) there follows. Continuity of tp v (9) and <pi n (0) over 
the sets ij+ n {ui G N*} and ij + n {«i N*} follows in an identical manner. 

Recall, now, that any isometry (p of the hyperbolic space extends to a 
smooth diffeomorphism of the conformal boundary of (J^,h), say %. This 
shows that on the set {v(9) > 0} the solution x := (f)ox will have a (/^-position 
function 

<Poo(0) = X(<A»(0)) ■ 

The fact that isometries extend smoothly to the conformal boundary further 
shows that the asymptotic behavior (4.5)-(4.6) is preserved under the action 
of isometries of the hyperbolic space, so that the map (p, (p) will satisfy the 
hatted version of (4.5)-(4.6) at points with positive v\. 

Let 9 be any point such that x( ( Poo(^)) 7^ tt mod 2tt. A straightforward 
analysis of (3.12)-(3.13) shows that one will recover (3.5)-(3.6) for the map 
(P,Q) at 9, with ^(9) = Vl (0) > 0. 

Set 

h+ ■= {x(<Poo{0)) * mod2vr}M 

V v ' 

{3 interval J around 9 such that \ ° ^ooIj = tt mod 2-7r} C ij+ . 

V „ ' 

U+- 

Then /j + is clearly open in Suppose that 6* G ij+ is such that there are 
no points of in a neighborhood of 9, then is constant and equal 
to 7r mod 2-7T on that neighborhood, hence 9 G /*+-. It follows that Ij+ is 
dense in 

On the function — 7r avoids integer multiples of 2ir, and con- 
tinuity of v\{9) = v\{6) > on ij + follows. Similarly one obtains a new 
continuous position function Qoo by using (4.1) with (p^ there replaced by 

On the other hand, at points at which ip^ = tt mod 2ir and v\ g" N* we 
have from (3.12)-(3.13) and (4.5)-(4.6) 

P(t,6) = Vl (9) In |t|+poo(g) + In ^1 + feM^ + (1) , (4.7) 

QM) = -^)-^°°)^(0) + O (l). (4.8) 

Analogous equations hold with supplementary In \t\ terms for tij £ N*, com- 
pare (3.4). This shows that (3.5)-(3.6) hold again, with a negative v\. Pur- 
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ther v\ is continuous on It follows that ij + contains an open dense set 

on which v\ is continuous. 

Consider, finally, points at which v\(9) < 0. Equation (4.3) leads to 

P (t, 6) = (-P+ ^ + ^ +O^)) (t,0) 
= -|«i(0)| ln|t|-P 00 (^) + o(l) . 
Inserting this into (3.13) yields 

sin(p) = -^—Q = 2e 2P -(l + o(l))\t\ 2lvim Q , 
sinn p 

so that s\n{ip{t,9)) goes to zero as t does. Equation (3.12) shows that we 
must have tp(t, 9) — >t_»o 7r mod 27r and one obtains, again modulo 2ir, 

f |t|2M*)| 2e 2Pao(*) (Q 0O (6») + (1)) , Wl _N* ; 

¥>M)-tt+< w2\vm 2 ^~m(Q rn (9)]n\t\ + Q x (0) + o(lj) , ^i(#) € —AT* 



(For v\(9) G — Af* we have used the hypothesis that In |i| xo(l) remains o(l).) 
The calculations done so far show that (3.5)-(3.6) hold for all 9 G [a, b] for 
the map (pox. A repetition of the arguments given on ij+ justifies continuity 
on an open dense subset of and the proposition for <ft\ easily follows. 

The result for % follows immediately from the calculations above, using 
the fact that for maps in % all level sets of Qoo form a finite collection of 
points. □ 



5 A symmetric hyperbolic system 

Let us set 

Pt = h, Pe = 9i, 
e P Qt = h, e p Q e = 92- 

Equation (1.3) takes then the form of the following first order symmetric 
hyperbolic system: 



d t 



( h \ 




(° 





1 


\ 




( h \ 




/ 


fi- 9 2 2 -^ 


\ 


h 













1 


de 


h 


+ 




-hh + 9192 - t 








1 











91 









\ 92 J 




V 


1 





J 




\ 92 ) 




V 


fi92 ~ 91 j 2 


/ 



(5.1) 
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The new unknowns f a ,g a agree with the coefficient functions of 1— forms, 
■ip^, of Christodoulou and Tahvildar-Zadeh's work on spherically symmetric 
wave maps [3], when an appropriate trivialisation of the bundle of vectors 
tangent to the hyperbolic space has been chosen: Indeed, if we set 

ei = dp , e 2 = e' p d Q , (5.2) 

then, in view of (1.2), the e a 's form a globally defined /i-orthonormal frame 
with constant structure coefficients (and thus constant connection coeffi- 
cients), and 

fa = h(e a ,X t ), g a = h(e a ,X e ). 

We consider solutions defined on domains of dependence fi(a, b, to), de- 
fined in (1.5). By Proposition 1.4 we have 

sup (|t/i| + |t/ 2 | + |t0i| + M)<oo. (5.3) 
Ct(a,b,to) 

Throughout this work the value of various irrelevant constants may 
change from line to line. 

Since \Xq\ 2 = g\ + g 2 , Proposition 1.6 implies 

lim^o f b t t 2 (g 2 1 +g 2 2 )d9 = 0, (5.4) 

Ja+t 

b-t 

t(gl + g 2 2 )deeL\[t Q ,0]) . (5.5) 

-+t 

We further note that by (1.4) we have 

D^Xv = (d^P - e 2P d^Qd u Q)d P + {d„d v Q + Q^P U + Q v P^d Q 

= (d^P - e 2P d^Qd u Q) ei + e p {8 ll d v Q + Q^P U + Q u P^)e 2 , 

(5.6) 

so that 

D e X e = (d e gi - gl)ei + (d e g 2 + gig 2 )e 2 . (5.7) 
It then easily follows from (5.3)-(5.5) together with Proposition 1.6 that 

lim^o f ' t 4 {{degif + (d e g 2 ) 2 ) d6 = , (5.8) 

Ja+t 

" 1 1 {(d e9 i) 2 + (d e92 ) 2 ) d9 G L\[t , 0]) . (5.9) 

-+t 

It turns out that we also have 



/ 

J (l 



/ 

J (l 
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Proposition 5.1 



lim t 



rb-t 

,o / t 2 /* 

Ja+t 



rfjffl = o , 



-b-t 



Proof: Let 



/ tfldO £ L\[t Q ,0]) . 

Ja+t 

F(t)= f 'tfide, 

Ja+t 



(5.10) 
(5.11) 



then 



dF 
~dt 



rb-t 

t( 9l - f 1 )(t,b-t)-t{g 1 + f 1 )(t,a + t)+ / t(f. 

Ja+t 



I " 9l) 



Now, tfi is bounded, hence so if F, and by integration of the last equation 
we obtain 



-t r b-t 



f f tfi d9dt 

J to J a+t 



nt rb-t 

~ / / 

J to Ja+t 



tg\ dOdt + C 



for all to < t < 0. Equation (5.5) together with the monotone convergence 
theorem imply (5.11). In order to prove (5.10) we calculate 



dtj a+t ^ 



+ 



-t 2 fi(t,b-t)-t 2 fi(t,a + t) 

rb-t 

/ h {de92 - f\h + 5132) 

Ja+t 



b-t 



< C{1+ I {\t\ 3 (de92f + \t\(fi + gi+g 2 2 )) 

I a+t 



(5.12) 



The function on the right-hand side of the last line is in L 1 ([to,0]) by (5.5), 
(5.9) and (5.11). Integrating (5.12) one concludes that the limit at the 
left-hand side of (5.10) exists. If this limit were different from zero (5.11) 
couldn't hold, whence the result. □ 



We are ready to prove now: 



24 



Proposition 5.2 For k > we have 

lim^o f * \t\ 2{k+1) i(d k e gi? + (fe) 2 ) d9 = 0, (5.13) 

Ja+t V 7 

f * If fe+1 ((fe) 2 + {d k e92?) d6 G ^([io, 0]) , (5.14) 

Ja+t K 7 

lim t _ f 1 \t\ 2 ^\d k e f 2 ) 2 de = , (5.15) 

Ja+t 

f \t\ 2k +\d k e f 2 ) 2 d6 £ L\[t Q ,ti\) . (5.16) 

Ja+t 

Proof: The cases k = have already been established, as well as (5.13) and 
(5.14) with k = 1. A simple induction argument, using the formulae (1.4) 
for the Christoffel symbols, shows that 

D k e X e = (d k g 1 +F k (d k - 1 g,...,g)y i +(d k g 2 +G k (d k - 1 g,...,g)y 2 , (5.17) 

where g denotes both gi, g 2 , while the -Ffc(-)'s and Gjt(-)'s are polynomials 
in the variables <9™g, < m < k — 1 with the number of derivatives and 
factors in each of the terms d^g-d^g... d l g n g satisfying 

n 

n>2, 5^(i m + l) < k + 1 . 

m=l 

This, together with Proposition 1.6, proves (5.13) and (5.14). 
Next, for k > 1 we compute 

k 

d k e f 2 = Y,C{i,k){DlX t ,D k e - l e 2 ) (5.18) 

i=0 

since f 2 = (Xt, e 2 ). By induction we obtain 

D k e 2 = (d k ~ l g + F^d^g, ...,g))e 1 + (d^g + G'^d^g, ...,<?)) e 2 , 

(5.19) 

where the F^_ 1 (-)'s and G' fc _ 1 (-)'s have the same property as described above. 
Thus 

k 

\t\ k+1 \d k f 2 \ < ^C(i,fc)|t| i+1 |DjX t ||t| fc - i |D*- < e 2 | 

i=0 

k 

i=i 
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The proof is completed by combining Proposition 1.6 with Equations (5.13)- 
(5.14). □ 



6 Existence of a velocity function 

Numerical experiments (see [2] and references therein) suggest that the limit 

v(9) :=lim\tX t \(t,9) (6.1) 

exists, and is a continuous function of 9 except for a "small" exceptional set 
of 0's in S 1 . We set 

f^t-rcg := {9 € S 1 such that the limit (6.1) exists } , (6-2) 

so that v is a well-defined function on f2t_ reg . The existence of this limit 
is useful when analysing the geometry of the associated space-time. Now 
Proposition 5.1 shows that t 2 f 2 goes to zero in L 2 as t tends to zero, and 
since 

\tx t \ 2 = t 2 (f! + fi) , 

the whole information about the limit (6.1) is contained in tf±, except pos- 
sibly for a negligible set. 

6.1 Existence of a weak velocity function v wea k 

We have the following: 

PROPOSITION 6.1 There exists f wca k G L°°(S' 1 ) such that for any p G (l,oo) 

\t\fl(t, ^ ^weak , 

where — 1 denotes weak convergence in L' P {S ). 

Remark 6.2 In the proof of Theorem 1.3 we establish existence of an open 
dense set Cl C S 1 such that t> wca k has a smooth representative v on Cl, with 
pointwise convergence to u on O. 

Proof: Let Vi{9) = 2~ % f\(— 2~ % , 9), then the sequence Vi is bounded in 
L 2 and therefore there exists f wea k and a subsequence Vi j which converges 
weakly to i> wea k- Let 4> be any smooth function on S 1 , we have 

d t [ tM = [ {td egi -t(f 2 -g 2 2 ))^ 
Js 1 Js 1 

= I {-tgM - t(f 2 - g 2 )4>) . 
Js 1 
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Integrating one finds 

I hh{t u o)<t>{p)dB - f t 2 h(t 2 ,e)<i>(e)d0 

Js 1 Js 1 

= I 1 I {-tgM - t(fi - gDcf) dOdt 

Jto Js 1 



Setting t\ 



t 2 J s 1 

and letting j go to infinity one obtains 



v weak (9)(f>(6)d6 - / \t 2 \fi(t 2 ,e)<f>(0)dB 
s 1 Js 1 



< f I (\tgid e <t>\ + + &)</>) dOdt . 
Jt 2 Js 1 

Since the integrand in the last line is in L 1 ([to,0] x S 1 ) we obtain 

hm/ \t\f l (t,e)<f ) (9)do= [ v weak {e)<t>{e)do , 

*->o Js 1 Js 1 

so that converges to f wea k hi the sense of distributions. Weak conver- 
gence in LP follows by elementary functional analysis, using the fact that 
smooth functions are dense in LP for p £ (l,oo), with p' - the Holder 
conjugate of p. Lower semi-continuity of the norm with respect to weak 
convergence implies 

ji) < lim ||2-^/ 1 (-2^,-)IIlp(51) < 2tt sup \th\ , 



\\v wea k\\Lp(Si) S Hm p " j Ji{-^,-)\\lp(si) 

[to,0)x5i 



so that 

||v W eak||L°°(S 1 ) < 7T SU P ll^wcakllLPfS 1 ) < SU P ■ 
27r pg[2,oo) [to,0)x5i 

□ 

The information contained in i> wea k seems to be very poor. For instance, 
since i> we ak is defined only almost everywhere, one can imagine situations in 
which f wea k has a smooth representative, but nevertheless the dynamics has 
very rough features at some points. This actually happens in the solutions 
with "spikes" discussed in Section 3. In those last examples one has point- 
wise convergence of \tXf \ everywhere; this suggests that even an exhaustive 
understanding of the properties of the velocity function might not be enough 
to understand the dynamics of the Gowdy models. 
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6.2 From space-time integrals to pointwise velocity 

We introduce 

Sgf (t)=g 1 (t,ij;-t)-g 1 {t,ijj + t) , (6.3) 

Ofl-reg := {V> e 5* 1 such that Jgf G L 1 ^, 0)) and 

/ (/| + <7?+5l)^^<ooi • 
^ t ° (V>) J 

Our aim in this section is to present an integral criterion for pointwise ex- 
istence of a velocity function v on Qe- rcg ; this will be used later in this 
work. 

Theorem 6.3 For every ip £ Qg_ reg there exists a number v(ip) £ M. such 
that 

lim c o^\tX t \ = v(ip) , with lim c o ^\tf 2 \ = (6.5) 
(recall that lim^o ^ has been defined in (2.4) ). This implies in particular 

^0— reg ^— ^t— reg • 

Moreover for ip £ f^-reg we /ictve 

lim c? ()W l t2 A^e| = lim^^l^DeXtl = lim^ \ tX e \ = . (6.6) 



Remark 6.4 We remark that for those points ip £ Oe_ reg for which v(ip) ^ 1 
we have curvature blow-up in the associated space-time. 

Proof: We have 

= d e g 1 -(l + 2tf 1 )(f 2 + g 2 ) + 2tf 2 d 6 g 2 -2tg 2 d e f 2 

+4tf 2 g ig2 , (6.7) 
a {tf 2 ) = d e g 2 + hf 2 + g 1 g 2 -2tf 2 d e gi + 2tg 2 defi 

+tf2g 2 2 -tf 2 g 2 1 -tf! + tfff 2 . (6.8) 

In particular 

□ (t/ 2 ) 2 = 2tf 2 o(tf 2 ) + 2(d t (tf 2 )) 2 -2t 2 (def 2 ) 2 

= 2tf 2 n(tf 2 )-2t 2 (d e f 2 ) 2 + 2t 2 (d e g 2 + fl-g 2 2 ) 2 (6.9) 
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using (5.1). It follows from (2.2), (6.10), Propositions 2.1 and [4, point 
b) of Lemma 3.4.1] that for V & ^e-reg the right-hand side of (6.9) is in 
L 1 (C'£(V ; )). The dominated convergence theorem applied to the usual inte- 
gral representation of solutions of the one-dimensional wave equation, 



1 rt px+t—s 

u(t, x) = u(t, x) + - / / a u (s, 9) dOds , 

2 J s =t Je=x-t+s 



where u is the solution of the free wave equation with the same initial data, 
shows that the limit lim^o ^(tf2) 2 exists. This limit has to be zero, oth- 
erwise the integral condition on /2 in (6.4) wouldn't hold. If Sgf is in 
L 1 ([to,0)), then the right-hand side of (6.7) is also in L 1 (C t ° (V')), and by a 
similar argument lim^o (^)tfi exists. □ 

We close this section with the following remark: 
Lemma 6.5 In the definition of Qg^ rcg , the condition 

/ f 2 2 dtd9< oo can be replaced by / \def\\dtd9 < oo . (6.10) 
Similarly, the condition 

/ \5gf\ dt < oo can be replaced by / \d$gi\ dtdO < oo . (6.11) 

J to ' JC° W 



Proof: Let 



then 



F(t)= f 1 h(t,6)d6, 



^ = -f 1 (t^-t)-f l {t^ + t) + f^^ + ^ + fl-g 2 ^ 

= -f 1 ( t ^-t)-f 1 (t^ + t) + dgi(t) + J^J(^ + fl-^j . 

Integration by parts gives the identity 

1 /"^~* ib - 9 

-r u{9)d9 = u{6-t) + u{9 + t)+ --—u'{9)d9, (6.12) 

1*1 Jib+t Jib+t \t\ 
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so that 

As F is bounded, integrating in t gives (6.10). Equation (6.11) is obvious. 
□ 



7 Power law in Sobolev spaces 

As discussed in the Introduction, and as will be proved below in detail in any 
case, a power law inequality (2.1) implies existence of a velocity function. 
It turns out that one strategy for establishing (2.1) is to derive a power-law 
for t-weighted Sobolev norms. This is done in this section. 
It is useful to introduce the following quantities 

Hi= sup |£/i| , [i 2 = sup \tf 2 \ , A 2 = sup \tg 2 \ ■ (7.1) 

Q(a,b,to) Q(a,b,to) Q(a,b,to) 

These are finite by (5.3). (Recall that Cl(a, b, to) has been defined in (1.5).) 
Let to < t < 0, a < b. Define the k— th order energy E^it) by 



E k (t) = \A 2k+2 £ ({d k eU? + fe) 2 ) dO. (7.2) 

J a+t o 



r b-t 

la+t i=l,2 

Proposition 7.1 1. If 



/?:= sup 1^1 + — + — < 1 , (7.3) 
n(o,6,to) \ 2 V 1 + 4 I^2| 2 + 1/ 

i/ien 

#k(t) < C(t ,k,a)\t\ 2a , fc>l, (7.4) 

where 

a = l-/3. (7.5) 

2. Similarly if 

$:= sup max(l-|t|/i,|*|/i)<l (7.6) 

fl(a,b,t ) 

then (7.4) ZioMs mi/t a = 1 — (3 if [i± < 1/2, and a - any number strictly 
smaller than 1 — fii if [i\ > 1/2. 
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Remark 7.2 Recall that, at fixed a and b, the constants fii and Aj depend 
upon to- I n several situations of interest (42 and A2 will tend to zero as to 
tends to zero, in which case the essential restriction in (7.3) is that fii be 
smaller than one sufficiently close to the singular boundary t = 0. 

Remark 7.3 We note that (7.3) will hold under the slightly stronger but 
simpler condition 

sup (| t/l | + Ml + M) <1 . (7 . 7) 

n(a,b,t ) V 1 V2 / 

Proposition 1.4 or Remark 1.5 can be used to replace (7.7) by a condition 
on initial data using the Cauchy-Schwarz inequality 



< 1 

t=t 



< y/7^\tX t \* + \tX e \ 2 
PROOF: Differentiating E^it) in t, using the field equations and integrating 
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by parts one has 



dE k (t) 
dt 



t 2k+2 £ ( {d k {f ._ gi)) 2 {tjb _ t) + {d k {f . + gi)) 2^ a + t ^ 



=1,2 

rb-t 



- {2k + 2) f \t\ 2k+1 ((defi? + fe) 2 ) 

Ja + t i=l,2 

+ 2|tp+ 2 T ' • - gl - f -}) + d k g 2 ■ d k e (f l92 - 9l f 2 ) 

Ja+t v T 



+ d$f 2 -d k e (-f 1 f 2 +g 1 g 2 -^)) 



d6 



<-2k f~* \t\ 2k ^ Y, (^/i) 2 " (2fc + 2) \t\ 2k ^ £ ( d 09i) 2 

Ja+t i=12 Ja+t i=12 

(7.8) 

/ (hd k ehd k e h ~ 92d k e hd k g 2 - f 2 d k g 2 d k 9l (7.9) 

Ja+t K 

+ 92d k e9l d k f 2 - h{d k f 2 f + h{d k g 2 f) 

(7.10) 

+ \t\ 2k+2 f 1 d k u. y, c(ij,k)diu-di u . (7.H) 

J a+t i i a—u 



b-t 

+ 2 | t |2fe+2 

I a+t 



rb-t 



'«+' i+j=k 

i,j>0 



In (7.9) and (7.10) we have collected all those terms which contain undif- 
ferentiated functions fi or g^. In (7.11) we denote (fi,gi) by u and the 
C(i,j, fc)'s are the coefficients of A:— th binomial expansions; we will ignore 
those coefficients and replace them by an overall constant from now on - 
this is sufficient for estimation purposes. Note that when k = 1, then (7.11) 
does not appear. Mixed terms of the form daf\daf 2 and d^g\d^g 2 are es- 
timated in the obvious way using 2ab < a 2 + b 2 . To take advantage of the 
different factors in front of the integrals appearing in (7.8) the mixed terms 
dgfidggj are estimated using 2ab < 2- + ab 2 . Absorbing all the terms from 
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(7.9) and (7.10) into those appearing in (7.8) we obtain 



^ < - \t? k+1 f (2k- 2|i/i| - \t f2 \ _?M) {d k e de 

M Ja+t * i=12 

- \t\ 2k+1 f b \2k + 2- 2\th\ - \tf 2 \ - 2a\tg 2 \) £ {d k e gtf 00 

Ja+t i=l,2 
b-t 



_|_ U|2fc+2 



C(k) f d k e u- d e u ■ % u ■ 

J a+t „• i _•_>. 



i,j>0 



For k = 1 this reads 



^ < - \t\ S [^(2 - 2\tfi\ - \tf 2 \ - 2 -M) J2 (d k d0 
dt J a+t i% 

-|t| 3 I {A-2\tf 1 \-\tf 2 \-2a\tg 2 \)Y,{dem?d6 . 

J a+t ■ i o 



It should be clear from what follows that the choice of a = a(t, 9) which is 
optimal for our purposes is that of equal factors in front of the sums, namely 
a = (-y/l + 4|t<72| 2 + l)/(2\tg 2 \). Choosing this value of a leads to 

dE x (t) 2a 
~dT ~ ~ |t| El ' 

with a as in (7.5). This shows that d((—t)~ 2a E\) / dt < 0, and by integration 
one obtains 

' t 



Ei(t) < 



to 



2a 



Ei(t ). 



(This inequality holds whatever the sign of a, but for a < it does not 
carry any new information.) 

To cover the case (7.6) we shall need a Lemma: 

Lemma 7.4 Let 

rb-t 



F(t)=f Iff + si) 

Ja+t 



I a+t 

Under (7.6) we have 

t ~' 2 '' 
F(t) < F(to) - 
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Proof: We have 

dF(t) 
dt 



< 



~{h ~ 92?{t, b-t)-{f 2 + g 2 )\t, a + t) 

l r l Ja+t 

\t\ ' 



and the result follows by integration as before. 



(7.12) 



□ 



Returning to the estimation of dEi/dt, assume that (7.6) holds and 
consider any of the terms of the form f^doudgu or g2dgudgu in (7.9) and 
(7.10); they are estimated as 



b-t 



a+t 



\t\ A f2d e ud e u 



< 



< 



< 



b-t 

a+t 
b-t 



\ -Mt\d e u\ 

Ja+t 4 



2 + ^1*1/1 



I a+t 

eE 1 (t) 
4\t\ 



C'(e)\t 



1-2/3 



where in the last line we have used Lemma 7.4; similarly for g 2 . It follows 
that 

dEi(t) 2(l-^i-e) 



dt 



< 



-^ 1 + C'(e)|t| 1 -^, 



Multiplying by l^i) 2 ^ w e ~> and integrating in t, one obtains 



Ei(t) < C(a) 



2 a 



with a = n\ if fx± < 1/2, and a - any number strictly smaller than 1 — fii if 
Mi > 1/2- 

The cases k > 2 are established by induction: Suppose, thus, that (7.4) 
holds for k = m — 1; we have already shown that it holds for k = 1. Then 
the terms from line (7.11) in ^^.^ are estimated as 



1*1 



2m+2 



rb-t 



' a+t 



d r e n u- diu-diu < C(t )\t 



m+l 



i+j = m 
i,j>0 



b-t 



u 



a +* i+j — m 

0<i,j<m 
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letting C(to) = sup^g) |i| m+1 |<9™ii|, which is finite by (5.3). Using the in- 
duction hypothesis for i,j<m on d l u ■ d^u, we get 

r b-t rb-t 

\ ^ I i 1 ?" -J— 1 I r\i I I i 11' —I— 1 I r\ ( ! 

3 



\t\ m+l [ ' e ifljui • i^ui = ir 1 / " E itrw-w'+w 

J a-\-t i+j= m J a-\-t i+j= m 

0<i,j <m 0<i,j <m 

<\t\~ x f'* E |t| 2i+2 (^n) 2 + |t| 2 ^ +2 (a>) 



2a-l 



< C(m)\t\ 
It follows that 

^ < - 2(m + °- 1) K > ( t ) + C7fa,m)| t |^ . (7.13) 

Multiplying by |^.| 2 ( m + Q - 1 ) on both sides and integrating over (to,*) in t, 
we obtain 

EUt)<\U 2{m+a ' l) E k {t ) + C{t Q ,m)\t\^ m+a -^ f\s\ l - 2m ds 
< C(t ,m)|t| 2a , 

as claimed. □ 



Corollary 7.5 Under the conditions of Proposition 7.1, there exists a con- 
stant C such that 



f \ 2 \Xg\ 2 de = f \ 2 (( gi ) 2 + (g 2 ) 2 ) d9 < C\t\ 2a . (7.14) 

J a+t J a+t 

If 



(71 := inf \t\h > , (7.15) 

fl(a,b,to) 



then we also have 



b-t 

t 2 f$d6 < C\t\ 2a ' , (7.16) 

a+t 



with a! = a if a < o\, or a' - any number smaller than o\ otherwise. 
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Proof: We calculate 
d ^ 



dt J a +t 



[ t 2 g 2 = - t 2 g 2 (t,b-t)-t 2 g 2 (t,a + t) 

Ja+t 

+2t f ' g 2 + 2 f t t 2 g2(d e h + f 1 g 2 -gif: 

Ja+t Ja+t V^v^ ^v— ' ^-v^ 

I II III 

9 rb—t rb—t -i rb—t 

< -7-7 / t 2 g 2 + e / \t\g 2 + - / \t\\d e h) 

\ l \ Ja+t Ja+t e Ja+t 



+t/:>f£/(^ 



II III 

< 2 - 2/ii - e'/z 2 ~ e /*'' ' ,-. , , ,m,i2«-i 



/ t 2 52 2 + C(e;e ')| t p 



1*1 ^a+i 

(7.17) 

Choosing e, e' appropriately and arguing as in the paragraph following (7.13) 
one obtains the bound for * g\. A similar, but simpler, calculation with 



c/2 replaced by g\ proves (7.16). In order to establish (7.16) we note that 
r b-t 

la+t 



if #8 = -t 2 f!(t,b-t)-t 2 f 2 (t,a + t) 

at Ja+t 

+2 J h \ 2 f2(dm-f^+9m) 
a+t I II in 

< e f %|/f + - f t \t\ 3 (d e g 2 ) 2 

Ja+t e Ja+t 



I 

\ l \ Ja+t \ l \ Ja+t 



< 



II III 

2(n - e '- b - 1 



V /" t2 f' + C(e)\t\ 

\ T \ Ja+t 



2a- 1 



and we conclude as before. □ 
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8 Pointwise power law 

We start with the following observation: 

Lemma 8.1 (i) Suppose that there exists a > such that 



rb-t 

/ t 2 \X e \ 2 aie < C\t\ 2a 

Ja+t 



I a+t 

Then for all multi-indices v we have 

l ' b ^ t /■,,.„ \ 2 

I a+t 



f b 1 (\t\\ v \ +1 \D v X e \) 2 dd < C{u)\t\ 2a , (8.2) 

Ja+t ^ ' 



(ii) If\tX e \ < C\t\ a p then we also have 

sn Vn(aM) \t\\^\D^X e \<C'(u)\tr P . (8.3) 

(hi) If the constant a in (8.1) satisfies a > 1/2 then (8.3) holds with a p = 
a -1/2. 

Proof: Equation (8.2) follows by a straightforward adaptation of the proof 
of [4, Proposition 3.3.1]. Equation (8.3) is obtained from [4, Remark, p. 73]. 
To establish point (hi) it remains to show that 

sup \t\\X \ < C\t\ a - 1 ' 2 . (8.4) 

Cl(a,b,to) 

That last inequality is obtained by applying the interpolation inequality 
no ii , r<\\f& ii 3 / 4 ii n 1 / 4 

\\O0U\\L°°([a+t,b-t]) S ^\\ e u \\ L ^[a+t,b-t])\\ U \\L^([a+t,b-t]) 

(see, e.g., [1, p. 94]; the condition there that u vanishes on the boundary is 
not necessary) to the functions tgi, i = 1, 2. □ 



Thus, power-law blow-up in Sobolev spaces implies a pointwise one if the 
decay rate is larger than 1/2. The unpleasant feature of the above argument 
is the loss of 1/2 decay rate in point (hi) of Lemma 8.1. This can be avoided 
by working directly with L°° norms, as follows: Consider two helds f,g 
satisfying the symmetric hyperbolic set of equations 

dtf-deg = S f , 

dtg-def = S g , (8.5) 
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set 

Ttt[f,g] = T ee [f,g] := \(f 2 + g 2 ) , T te [f,g] = T et [f,g] := fg , (8.6) 
and define 

Mf,g] = d v {T l / , \f,g]). (8.7) 

Writing T^ u for T^ u [f,g], etc, we have the identity [4, Equation (3.2.5)] 

r 



T tt (ti,ei) = ~ jT fo t + jo)(t,e 1 + t 1 -t) + (j t -jo)(t,o 1 + t 1 + tj) 

+- ((T tt + Tto)(t, 0i + ti - t ) + (T tt + T w ){t, 1 + t 1 + t ) 



dt 



For (8.5) the j-terms appearing in (8.8) read 

-\(k + je) = \{f + g)(s f + s g ) , 
-\{k-je) = \{f - g)(s f - s g ) . 

If we let (f,g) = we obtain 



(8.8) 

(8.9) 
(8.10) 



~Ut+je) = ^d k f 1 + d k gi )(^ + d k (/ 2 -52)(/ 2 +52) 



1*1 



\(jt-je) = \{d k eh-d k egi ){^ + d k e 



{h - 52) (/2 +52) 



(8.11) 



(8.12) 



Similarly for (/,#) = (d%f 2 ,d%g 2 ) one has 



\t\ 



&e (h + 9i)(f2-92) 



(8.13) 



--(.;/ -./V) - " 



1*1 



(/l-5l)(/2+52) 



(8.14) 



We define 



sup {T tt [dlh,d$ gi ] + T«[^/2, fe])(s, 9) . 

t <s<t, ip+s-t<e<i>-s+t 

(8.15) 
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It is useful to note that 



Define 



Tu[f,g] = \((f-g? + {f + g?) 



H± = sup \t(fi ±gi)\ ■ 

il(a,b,to) 



(8.16) 
(8.17) 



We have the following pointwise equivalent of Proposition 7.1: 
Proposition 8.2 Suppose that 

sup L- + \tfi + t 9l \ + + \th + t 9l \) 2 + 4\tf 2 + tg 2 \ 2 ) < 4 , 

Q(a,b,t ) v 7 

(8.18) 

sup (n + + \th - t gi \ + vW + \tfi ~ tgi\) 2 + 4|t/ 2 - t<72| 2 ) < 4 . 

n(o,6,t ) v 7 

(8.19) 

T/ien £/iere exisi constants C, a p > so £/ia£ we /ictve £/ie inequality 



\tX e \ < C\t\ 



(8.20) 



Proof: We start by deriving an integral inequality for E\ using (8.8). Let 
S denote the sum of (8.11) and (8.13) with k = 1: 

S = -(deh+d e9l )\ 



+ d e (f 2 -92){h + 92, 



. 1 ,r, , . a , ( d e f 2 + d e g 2 + d e f 2 - d e g 2 
+- 2 {def2 + d e g 2 ) ^ ^ 



1 

W\ 



(h-g2)(h+gl) 

{deh + degif + (d e f 2 + d e g 2 f + (d e f± + d e gi)(d e h - d e gi) 



+ (deh + d e g 2 )(d e f 2 - d e g 2 )^ + -(f 2 + g 2 )(d e fi + d e gi){d e f 2 - d e g 2 ) 



1 



(h + 9i) (deh + d 9 g 2 )(d e f 2 - d e g 2 ). 



(8.21) 



A formula for the sum of (8.12) and (8.14) can be obtained by changing g a 
to — g a in (8.21). We apply Young's inequality ah < |a 2 + ^-b 2 to estimate 
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the mixed terms (dgf + dgg){def — dgg), using a c which might possibly 
depend upon t and 9. With a little work one finds that the integrand in the 
first line of (8.8) can be estimated by the sup, over 6 in the relevant range, 
of the quantity 



1 

2\t\ 



1 . c l \tf 2 + tg 2 \ \ . . 2 
1 H (defi + degi) 



+ ( i + ^ + \ (\th + 1 91 \ + M±M) , W2 _ a 052) ^ 

+ {l+ C2\th-t 92 \ \ {deh _ dm? 



+ [l + ^ + U\th-tg 1 \ + ^^))(d e f 2 + d e92 ) 



(8.22) 



If all the factors in front of the derivative squared terms are strictly smaller 
than 2, say smaller than or equal to 2 — 2a p , then (8.22) is smaller than 
or equal to (4 — 4a p )-E'i/|t|. From (8.8) applied to E\ and from GronwalFs 
Lemma (c/., e.g., [4, Lemma 3.2.3]) one obtains an integral inequality for 
Ei, which translates into the inequality 

\t 2 def a \ + \t 2 d e g a \ < C\t\ a " . (8.23) 

It should be clear that an optimal estimate will be obtained in (8.22) if 
ci = ci (t, 9) is chosen so that 

C1I/2 + 52I = n- + 1 A + 01 1 + 



Cl 

similarly for c 2 . This gives 

ci\f 2 +g 2 \ < + + VG"- + l/i+Si|) 2 +4|/ 2 + 52 | 2 ) , 

C2I/2 - 52I < ^ (m+ + l/i - 51 1 + VG"+ + l/i-5i|) 2 +4|/ 2 -5 2 | 2 ) , 

with equalities if (f 2 + g 2 )(f 2 — g 2 ) 7^ 0, and the condition mentioned above 
will be satisfied if (8.18)-(8.19) holds. 
Equation (5.1) gives 

d t gi = 0(\t\ a ^ 2 ) , (8.24) 
and by integration along rays 9 = tp + Xt, A G [— 1, 1] one obtains on C^ (tp) 

\tgi\ < C\t\ a » . 
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Returning to (5.1) one finds 

d t 92 ~ fm = 0(\t\ a »- 2 ) , (8.25) 
which can be integrated to give 

> t 2 > *i > *o 92M) = efi fl{s ^ ds g 2 (t 1 ,iP) 

+ f * e-O h^) ds O{\u\ a v- 2 )du . 

(8.26) 

Since |/i| < (1 — 7) |ln |i|| + C one easily concludes that 

g 2 (t^) = 0(\t\ a ^ 1 ). 

Integrating in 9 from ip to 0, at fixed t, the desired estimate for g 2 on C° o (-0) 
is obtained using (8.23). Since tp was arbitrary in [a, b], and since all the 
constants were uniform in ip, the result follows. □ 



The main result of this section is the following: 
Theorem 8.3 Suppose that either 

sup ( m+ m + y )<i, ( 8.27) 

n(o,6,to) \ 2 V 1 + 4 I^2| 2 + 1/ 2 

or i/ta^ (8.18)-(8.19) /toW. Then the velocity u wea k a continuous repre- 
sentative v < 1 on [a,b], and the weak convergence in Proposition 6.1 can be 
replaced by convergence in sup norm to v. In other words, lim t ^ \tX t \(t,9) 
exists and is a continuous function on [a, b] . Moreover the solution satisfies 
a power law blow-up, Equation (8.20). 



Remark 8.4 Further information concerning the properties of the solutions 
considered in Theorem 8.3 can be found in Theorem 11.1 below. 



Proof: Under (8.27) Corollary 7.5 applies, so that the conclusion of Lemma 8.1 
point (iii) holds. It now follows from Proposition 8.2 that both under (8.27), 
or under (8.18)-(8.19), point (ii) of Lemma 8.1 applies, and thus there exists 
e < 1 such that \tDgXff\ + \tDgX t \ < C\t\~ e . Theorem 8.3 is now a straight- 
forward consequence of the following: □ 
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Lemma 8.5 Suppose that there exist positive constants C and a p such that 
on £l(a,b,to) we have 

\t 2 D e X e \ + \t 2 D e X t \ < C\t\ a ? . (8.28) 
Then there exists a continuous function v such that on Q(a, b, to) it holds 

\tX t \ 2 {t,9)-v 2 {9) <C'\t\ a ". (8.29) 

Further, for every if) such that (8.28) holds on C® (ip) we also have (8.29) 
onCHi,). 

Proof: We have 

dt\tX t \ 2 = 2h(tX t ,tDgXg) . (8.30) 
By integration we obtain 

\tX t \ 2 {t^9) = \tX t \ 2 (t 2 ,9) + f 1 d t \tX t \ 2 (s,e)ds = f 1 0{s a ^ 1 )ds . (8.31) 

It easily follows from this equation that 

v{6) := Km\tX t \(t,0) 

exists. By passing to the limit t\ — > in (8.31) we obtain on [to , 0) x [a, b] 

r0 



\v 2 



{9) - \tX t \ 2 (t,9)\ < J Cs^^ds < C'\t\ a p . (8.32) 



This shows that \tXt\(t,-) converges uniformly to v, and establishes conti- 
nuity thereof. The same argument applies on C® (ip) by integrating along 
rays 9 = ip + \t, A G [— 1, 1]; one easily checks, using (8.28), that the number 
v(ip) is A-independent. This, together with the result already established 
on [to,0) x [a, b], establishes (8.29) on fi(a, b, to). □ 

In terms of the (/, g) variables the Gowdy-to-Ernst transformation (3.20) 
takes a remarkably simple form: 

A = —fi - ~ > 9i = ~9i , 

A = ~92 , 92 = ~h ■ ( 8 - 33 ) 

This transformation immediately leads to the following counterpart of The- 
orem 8.3: 
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Theorem 8.6 Suppose that either 

sup (| tA + 1 | + M + (8.34) 

or 

sup (A_ + + 1 + t gi \ + + + 1 + ^i|) 2 + 4|i/ 2 + ^ 2 | 2 ) < 4 , 

n(a,6,t ) V 7 

(8.35) 

sup f A+ + \th + i - t9i\ + VG"+ + + 1 " ^i|) 2 + 4|t/ 2 - ^ 2 | 2 ) < 4 , 

f2(a,Mo) V 7 

(8.36) 

where 

A± = su p + 1 ± tffil • 

Q(a,6,to) 

T/ien i/iere exisfe a continuous function v such that 

\th\ 2 + \tg 2 \ 2 ^ v 2 , (8.37) 
uniformly in 9 . Further there exist constants C, e > such that we have 

\tf2\ 2 + \tgi\ 2 <C\t\ 2 ' . (8.38) 

In Section 12 below we will see how to iterate the Gowdy-to-Ernst trans- 
formation to obtain information on more general solutions. 

For further purposes it is convenient to restate the conclusions of Lemma 8.1 
as higher derivative estimates for the g a J s and / a 's: 

Lemma 8.7 (i) There exists a constant C such that 

\t\ k+1 \d k e f a \ + \t\ k+1 \d k g a \ + \t\ k+l \d k f a \ + \t\ k+1 \d k g a \ < C . (8.39) 



f * (\t\ l+k+1 d?d l e u) 2 d6 < C(l,k)\t\ 2a , (8.40) 

Ja+t ^ ' 



(ii) Under the conditions of point (i) of Lemma 8.1 we have 
r6 ~* / . . , . , , , \ 2 

la+t 

for every I, k > when u = gi, g 2 or I > 1, k > when u = f±, f 2 . 
(hi) If\tXe\ < C\t\ ap then we also have 

^l+k+i^QkQl^ < C (l,k)\t\ a " , (8.41) 

with the ranges of the (I, k) 's as in (ii). 
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Proof: For the purposes of the proof let u stand for any of g a , f 2 . 

(i) : From (5.17) we have 

{d k e gi? + (d k g 2 f < \D k e Xe\ 2 + Fl + G\ 

then |i| fc+1 |<9|g| < C inductively using part (ii) of Proposition 1.4. For 
|^|fc+i|^fej a | we com p U te by induction 

D k t X e = (d k 9l + F k (d k - l u, u)) ei + (d k g 2 + G k (d k ~\ u))e 2 , 

(8.42) 

D k t X t = (dtf 1 + F k (d k - 1 f,...,f))e 1 + (dZf2 + G k (d k - 1 f,...,f))e 2 , 

(8.43) 

where F k , G k , F k , G k have the same properties as described in the proof of 
Proposition 5.2 for F k and G k . The remaining inequalities follow as before. 

(ii) : The proof is identical to that of Proposition 5.2. 

(iii) : We consider the g a 's first. Letting e stand for the basis vectors e a 's 
of (5.2), it is sufficient to show 

sup \t\ i+j \D\Dle\ < C, (8.44) 

Q(a,&,to) 

then the assertion follows from 

I 

d k d l 6 g = d k d l 6 (X s ,e) =d k Y,C(l,j)(DlX e ,D l -i e ) 

3=0 

k I 

= E E C & k > 3){DiDix e ,D*- i D l - j e) 

i=Q j=l 

together with Lemma 8.1. Now 

F>e&\ = 92e 2 , D e e 2 = -g 2 e,\ , D t e x = f 2 e 2 , D t e 2 = -f 2 e x , (8.45) 
which implies that the D^e a , F)\e a are of the form 

D\e x = F fc _ iei + (d k t - l f 2 + Gfc-x) e 2 , D k t e 2 = (-d*' 1 f 2 + e, + G' k _ x e 2 , 

D k ei = F k _ iei + (d k - 1 g 2 + G fc _i) e 2 , D k e 2 = {-d k ' l g 2 + i^-i) ei + G' k _ x e 2 , 

with the Fk—is, G k -\s, etc., of a similar structure as in (5.17). By (8.39) 
and by the above expression we see that |i| fc |Dge|, |t| fc |D^e| are bounded. 
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For mixed derivatives of g, we write 

D k t D l e X e = D k (id l e9l + F,)ei + (d l e g 2 + C?,)e 2 ) 
fc 

= £ aj(^ + F,)A fc_i ei + aj(fe + Gi)A* 

i=0 

= (#^ + #F0ei + (#fe + #G,)e 2 
fc-i 

+ J>|4<7i + diF^Dt'e! + (did l e g 2 + d\G l )D k - i e 2 ; 



i=0 

above, and in what follows, we ignore constants arising from binomial expan- 
sions. From above expression we get \t\ k+l+1 \d k d l e g\ < C inductively using 
(8.39) and boundedness of |t| l |£)|e|. Finally (8.44) follows from writing 

D k t D l e e = D k t ((a^- 1 ^ + F^e, + (ad^g + (Vi)e 2 )) 
k 

= Y,(°dld l f 1 92 + d\F l _ x )D k t - i e 1 + (ad^g + d i Gi-i)D k ~ i e 2 
i=0 

(8.46) 

and then using boundedness of \t\ l+ ^ +1 \d\d :, e g\, \t\ L \D\e\. 

Let us turn now to the / 2 's. According to (5.18) we write 

l 

d l e f 2 = J2(D J e X t ,D l d - j e 2 ), 

3=0 

so that 

Q k Ah = Y J ^ Di t D i x ^ D t~ iD V j ^), 

i=0 j=0 

= j^iDlXuD^D 1 ^) + jy£(E\& 9 X t ,n!<-*ltf> e 2 ) □ 

i=0 i=0 j=l 

when I > 1, k > 0. Multiplying by we have 

k k I 

\t\ k+l+1 \d k d l e f 2 \ <^|t| fe+ ^p t fe -^e 2 |+^^|t|^ +1 | J D^X t | (8.47) 

i=0 i=0 j=l 
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from Proposition 1.4 together with (8.44). Using the expression (8.46) we 
get 

k— i 

\ t \k+i-i\ D k-i D i e 2 \ < \t\ k+l - 1 (\dTd l i l g\ + {d^F^l + |D*- i - m c 

m=0 

< c itr +m {\d?$ e - i g \ + i^i-ii + i^r^-ii) 

m=0 

from the boundedness of |t| fc |-D^e|. The claim on f 2 follows now from Lemma 
8.1 together with the assertions on g. Finally, the result for fi follows from 
what has been proved so far together with the equation 

dt(tfi) = t{d e gi + fl-g 2 2 ) . 

9 Stability of the (|, |, -\) and (1, 0, 0) Kasner met- 
rics 

The following result extends the singularity stability theorem for the (|, |, — |) 
Kasner metrics established in [4], by raising 9 the stability threshold there 
by a factor 6 3/2 /2 

Theorem 9.1 Suppose that 

sup tl {\X t \ 2 + \X e \ 2 ) (to,0) < \. (9.1) 

6e[a-t ,b+to] 1 

Then the solution is of power-law type. Further, the curvature scalar Ra^&R 01 ^ 5 
blows up on every causal curve with endpoint on {0} x [a, b] x S* 1 x S 1 . In par- 
ticular the associated Gowdy space-time is inextendible across the boundary 
{0} x [a, b] x S 1 x S 1 . 

Proof: We wish to apply Proposition 8.2. Let ji± be defined by (8.17), 
using the Cauchy-Schwarz inequality and some rather obvious estimations 
we have 

M- + + tgi\ + Vin- + \tfi + tffi|) 2 + 4\tf 2 + tg 2 \ 2 ~ 

< \/2V2 (//_ + \th + t gi \Y + 4\tf 2 + tg 2 \* 

< 2>/2yJn 2 _ + \th + tgi\ 2 + \tf 2 + tg 2 \> . (9.2) 
9 A similar result has been recently established by Ringstrom [18]. 
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This shows that both conditions (8.18)-(8.19) will hold if 

sup (\tf 1 -tg 1 \ 2 + \tf 2 -tg 2 \ 2 ) < 1 , 

Cl(a,b,to) 

and if 

sup (\ t f 1+ t gi \ 2 + \tf 2 + tg 2 \ 2 ) < 1 . 

Q(a,b,to) 

It follows from (1.8) that the last inequalities will hold if 

sup \t\ {\fi-gi\ 2 + \fi+gi\ 2 + \f2-9 2 \ 2 + \f2 + 92\ 2 ) < 1 , 

t=to,ee[a-\to\,b+\to\] 

which is equivalent to (9.1). The result follows now from the arguments in 
the proofs of Theorem 3.5.1 and Proposition 3.5.2 in [4]. □ 



Remark 9.2 We note that the Sobolev decay estimates of Section 7 lead to 
a similar somewhat weaker statement, with 1/2 replaced by -y/3/19 in (9.1). 
This can be seen as follows: Equation (8.27) will hold under the slightly 
stronger but simpler condition 

sup ( m + m + m<\. (9.3) 
n(o,6,t ) V 1 V3 / * 

The Cauchy-Schwarz inequality, 

\th\ + ^ + 7f^y i + i + ^Vlt/il 2 + \tf2\ 2 + w , 

together with point (i) of Proposition 1.4 show that Theorem 8.3 applies, 
and one concludes as before. 

The solution 

/i = l/|t|, h =91 =92 = 0, 

of (5.1) corresponds to the flat Kasner metric. Theorem 8.6 similarly implies 
complete control of the behavior of the solution for all data in a neighbor- 
hood of those for the flat Kasner metric. In this case the geometric inter- 
pretation is more complicated, because of occurrence of horizons. A further 
discussion of the latter can be found in [6]. 
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10 Behavior of power-law solutions at v = and 

v = 1 



We consider solutions on Q(a,b,to) such that 

\tX e \ < C\t\ a * , (10.1) 

for some a p > 0. It follows from [4, Remark 7.3] and Lemma 8.5 that the 
velocity function v exists. The aim of this section is to study the behavior 
of such solutions at points, or intervals, on which v vanishes: 

Theorem 10.1 Suppose that (10.1) holds with some a p > and consider 
any point ip G [a, b] such that v(ip) = 0. Then: 

(i) The restriction x := x\ c o ^ of x to C^ Q (tp) can be extended to an 

map from M 2 to J%2 ■ 

(ii) // lim c o ^tdffPt = for all j G N, then for i, k £ N we have 

\^ w df +l d k d P = lim C QW d? +1 d$Q = . (10.2) 

(hi) Further, if v vanishes on an interval [9i,0 r ], then the restriction x := 
x \n(6i,6 r ,t ) °f x t° ^(^/j @r, to) can be extended by continuity to a smooth 
map from M 2 to J^, with (10.2) holding for all tp G [9i,0 r ]. 

Theorem 10.1 says, in essence, that x behaves on C^ (ip) as if it arose 

from an AVTD^' Q) map defined on R 2 . We emphasize, however, that the 
extensions mentioned above might fail to coincide with the original map x 
away from C t ° (i/>) (for x), or away from £l(6i, r ,to) (for x). Such a situation 
could arise when x has an infinite number of smaller and smaller spikes 
accumulating at a point at which v(ip) = 0. 

Equation (10.2) says, roughly speaking, that P and Q can be thought of 
as smooth functions of and t 2 (rather than t) . This result is relevant to the 
question of extendibility of the associated metric across Cauchy horizons. If 
lim^o ^tdgPt = for a finite number of j's, there will be a certain number 
of i's for which (10.2) will hold. 

Using the Gowdy-to-Ernst transformation we obtain immediately the 
following counterpart of Theorem 10.1 at v = 1: 
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Theorem 10.2 Suppose that there exist constants C, a p > such that 

|t0i| + |*/ 2 |<C|*| Q * (10.3) 

and consider any point ip G [a, b] such that vi(tp) = 1. Then the functions 
(P, Q) := (P + In \t\, Q)\c° <ib) can be extended to an AVTD^o'^ map from 
IR 2 to J%2- If li m c t ° (i/>)tdoPt = for all j G N, then for all i, k G N we have 

lim c o o(V)) 5f +1 d$(P + In \t\) = \im c o oW d? +1 d$Q = lim^ W tf +1 Q = . 

_ ( 10 - 4 ) 

Further, if v± = 1 on an interval [0i,0 r ], then the restriction (P,Q) := 
(P + In |£|, Q)|n(6»(,er.,to) can be extended to a smooth map from M 2 to J%2, 
with (10.4) holding for all tp G [9i,6 r ]. 

Remark 10.3 The vanishing of the last term in (10.4) for all k > is 
somewhat surprising. We emphasise that the power-law condition (10.1) in 
Theorem 10.1 is justified for small initial data by Theorem 8.3, and that the 
condition (10.3) is justified for initial data near the flat Kasner by Theo- 
rem 8.6. 

PROOF of Theorem 10.2: The Gowdy-to-Ernst transformed map x sat- 
isfies the hypotheses of Theorem 10.1, and therefore (10.2) holds for the 
associated functions P and Q. The claim about P in (10.4) is straightfor- 
ward. From (3.20) we have 

d t d k Q = -\t\d k e (e 2p d e Q) , d k e +l Q = -\t\d k e (e 2p d t Q) . (10.5) 

Integrating in t and using (10.2) one obtains (10.4). □ 



An iteration of isometries and Gowdy-to-Ernst transformations, as in 
the proof of Theorem 1.1 below, allows one to control the behavior of x near 
points (0,V>) at which v(ip) G Z, the details are left to the reader. 



Proof of Theorem 10.1: We start with a lemma: 

Lemma 10.4 Suppose that there exists < a p such that (10.1) holds. Then 
there exists a function v\{9), with \v\\ = v, and constants C a {9), a = 1,2, 
such that on C^ Q (9) we have 



\th\ < C 2 (9) { 



\t\ a ", 

i»M + \t\ a p, 



v{9) < 0; 

< «i(0) / a p ; 



{ (l + \\n\t\\)\t\^ e \ v 1 {9) = a. 



(10.6) 



v 
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' \t\ a % ^i(0)<O; 
\t\fi- Vl (9) < 0,(9)1 \t\^W + \t\^, 0< Vl (6)^a p ; (10.7) 
[ \t\ a *, v 1 (9) = a p . 

The constants C a (6) are uniformly bounded on compact intervals on which 
vi is strictly positive, uniformly bounded away from zero. 

Remark 10.5 The examples discussed in Section 3 show that v\ is not con- 
tinuous in general. Further, the constants C a are not uniformly bounded 
near points at which v\ has discontinuities involving a change of sign. Sim- 
ilarly we do not expect uniformity at points at which v\ crosses 0. 

Proof: As explained in the proof of Lemma 8.1, Equation (10.1) implies 

\t\\ v \+ 1 \D v X e \ < C(u)\t\ a p . (10.8) 

We can use Lemma 8.5 to conclude that there exists a continuous function 
v such that 

\v 2 {6) - \tX t \ 2 (t,9)\ < C\t\ a P . (10.9) 

Equation (10.8) with v = 9 and v = t shows that the ^-derivatives conditions 
of Lemma 6.5 are satisfied, and Theorem 6.3 implies that there exists a 
function v\ with = v such that, for all ip G [a, b], 

lim c o oW \t\h = Vl (^) , (10.10) 

together with 

lim c o oW |t|/ 2 = 0. (10.11) 

Equations (5.1) and Lemma 8.7 lead to 

Wl/i) = \t\f 2 + 0(\t\ a ^), (10.12) 
dt(\t\f 2 ) = -ltl/i/2 + Odtr- 1 ). (10.13) 

Suppose, first, that 9 is such that v(9) = 0. It follows from (10.9) that 
the right-hand sides of (10.12)-(10.13) are 0(\t\ ctp ~ 1 ). The same is true for 
do(\t\fa) by Lemma 8.7, and by integration in t along rays T\(t) = (t,9 + Xt), 
A G [-1,1], we find on C t ° o (0) 

v {9) = => |t/i| < C\t\ a " , \tf 2 \ < C\t\ a " . (10.14) 

In general, integrating (10.12) shows first that i/f G L l ([to,0), and then 



m- J\s\f: 



t\f 1 (t,6) = v 1 (9)- \s\fj(s,0)ds + O(\t\ a >) . (10.15) 
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Similarly, (10.13) can be integrated to give 

o>t 2 >h>to M/ 2 (i2,0) = e-^^^ihihi^e) 

+ [ 2 e-Pu h{sfi)ds ^ ap -i^ du 
Jti 

(10.16) 

If vi(9) > one finds first from (10.16) that \tf2\ < C\t\ e , with e equal to, 
say, mm(vi(6)/2, a p /4). Plugging this in (10.12) one obtains that f\ — v\ = 
0(\t\ e ). Returning to (10.16) one is then led to (10.6) at (t, 0) with vi{9) > 0. 
Inserting (10.6) into (10.15) we arrive at (10.7) at it, 9), again for vi(8) > 0. 
Suppose, finally, that v\(9) < 0, then we rewrite (10.16) as 

0>t 2 >t 1 >t M/2^,0) = e^ h{s ' e)ds \t 2 \f 2 (t 2 ,9) 

+ [*e I Z M >M da O(\u\ a *- 1 )du. 
Jti 

(10.17) 

One readily checks that the integrand is in L 1 ([to,0)), and passing to the 
limit t 2 —> one obtains 

0>h>t IhlMtuO) = f e^ h{sfi)ds O{\u\^~ l )du. 

Jti 

(10.18) 

Arguing as before one obtains the result at (t, 9). Finally, the result for 
(t, 9i) G Ct (9) is obtained from the one for if, 9) by integrating d${tf a ) in 
9 from (t, 9\) to (t, 0), using point (hi) of Lemma 8.7. □ 



We return to the proof of Theorem 10.1. To proceed further we need 
better control of the 9 derivatives. Let be defined by (8.15) and con- 
sider a point tp sucri that v(if>) = 0. From (10.1), (10.14), together with 
Equations (8.8) and (8.11)-(8.14) one obtains 

/•*! 9 + 0(\t\ a p) 

E 1 (t 1 ,i/>)<E 1 (t , 1 />)+ - v 1 '- Ekfa^dt; (10.19) 

Jto m 

we have also made use of (8.16). Gronwall's Lemma (c/., e.g., [4, Lemma 3.2.3]) 
gives, for any e > 0, decreasing \to\ if necessary, 



E 1 (t,^) < C 



2+e 



(10.20) 
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This allows us to rewrite (10.19) as 

ExM) < EiM) + £ (li^'VO +0(\t\ a *- 2 - e j) dt , (10.21) 

which, together with Gronwall's Lemma, implies (10.20) with e = 0. It 
follows that on {ip) we have 



\tdef a \ + \td e g a \ < C , 

and (5.1) gives 

\td t9l \<c, \td t92 \ < c(i + iti 20 *- 1 ) . 

Integrating in t over rays 6 = if) + At, A G [—1, 1], one obtains 
sn = 0(ln |«|) , Isal^CMtl + ltl 2 **- 1 ). 



(10.22) 
(10.23) 
(10.24) 



We have thus recovered (10.1) with the exponent a p there replaced by 2a p if 
2a p < 1, or by a number as close to one as desired otherwise. Iterating the 
argument leading from (10.23) to (10.24) a finite number of times if needed 
one thus has 

\tf a \ + \tg a \<C\t\ l -\ (10.25) 

with e as small as desired. 

We continue by induction: suppose that for all e > and for all < i < 
k — 1 there exists Cj(e) 



Mf a \ + \tdlg a \ < d{e)\t 
From (8.8), (8.11)-(8.14) one obtains 



(10.26) 



to 



< E k (t ,iP)+ [ 
Jt< 



E k {t^) + C\t\- 2 -*E\ /2 {t,$j) dt 

dt, (10.27) 



h /2 + 2e^ „ . C,.,_ 2 _ 2e 



-E k {t^) + -\t\ 



leading to (10.26) with i = k (and with e replaced by 2e.). 
The equation 



d$gi 

V $92 J 



( 1 \ 

1 

10 

\ o i o o / 



do 



( d k eh \ 

d k eh 

dfigi 
V d k e92 J 



+ 



( d$(fl 

d k e (-hh + 9i92-ff) 


de{fi92-gih 



(10.28) 



52 



gives 

\d t (tdif a ) | + \td t dig a \ < Q(e)|t|- e . (10.29) 

Integrating in t over rays = if) + At, A G [—1, 1], from (10.26) and (10.29) 
one obtains 

\dj>9a\ < Q(e)ir e • (10-30) 

Further, the same integration argument shows that there exist constants 
v a ,i(ip) such that 

||Wa-« O) iW0| < C^ltl 1 " 6 • (10.31) 

We note that v a fi(ip) = by (10.25). From (10.28) one obtains now the 
equation 

d t {di\t\h) = \t\di(fi)+o(\tr) 

£ 0) V2, j (^2,i- j (^ ^ + 0(|tp e ) , (10.32) 

which is compatible with boundedness of d l e {\t\f\) if and only if the coefficient 
of l/\t\ vanishes. Suppose that we know that V2,j vanishes for j = 1, . . . , m — 
1, then the condition of vanishing of the offending term in (10.32) with 
i = 2m gives f2,m — 0, and induction gives the vanishing of V2j for all j . It 
then follows from (10.31) that 

m 2 \ < C 4 (e)|tr • (10-33) 
Integrating the equations for dtd l e g a one arrives at 

8i 9l = -v Vl+l {i,)\n\t\+ gi , t {i,)+0(\t\ 1 -') , (10.34) 
92 = 9 2,oW + 0(\t\ 1 - e ), (10-35) 
492 = -v 1 ^)g 2fi Wln\t\+g2,iW + 0(\t\ 1 -') , (10.36) 

for some constants g a ,i(ip)- Similarly, as in (10.18), 

0>h>t |ti|$/ 2 (ti,^) = f e^^^O^-^du 

= 0(|ti| 2 " £ ). (10.37) 
Integrating in 9 from tp to 6 and using (10.37) one obtains, on C° (■)/>), 

d k e f2(t^)=0(\t\ 1 -*). (10.38) 



53 



Analogously to (10.15) one has 

\t\dlf 1 (t,^)-v 1 ^) = J" 0(\ S ln\s\\)ds = 0(\t 2 ln\t\\) . (10.39) 

An iterative repetition of the arguments given gives a full asymptotic ex- 
pansion of the solution on C° Q (ip). It should be clear that the expansions for 
the derivatives obtained in this way behave as though they arose from an 
AVTD^' Q) map, and the usual extension arguments can be used to provide 
an AVTD^ > ' ( ^- ) extension of the restriction x\ c o ^ of x. 

Suppose finally that all the vij(ip) vanish. We note that it follows from 
the argument above that this will be the case when v = on an interval 
containing tp. It is then straightforward to show that all the c^/ a 's and 
the dggaS can be extended by continuity to the tip (0,ip) of C^ip). By an 
abuse of notation we will denote the value at the tip of those extensions by 
<9 Q / a (O,-0), etc. Equation (10.37) gives, for all k, 

d k f 2 (0^)=0. (10.40) 

Equation (10.39) implies 

d e fe /i(0,V)=0. (10.41) 

Arguing similarly one obtains that any mixed t and 6 derivatives can also 
be extended by continuity to the tip (0,ip) of C° Q (ip). For example, when 
vijW = (10.39) can be rewritten as 



\t\d k e h{t^) 



e \U2\ + Q (|t| 3 ) , (10.42) 



which shows that dafi/t extends continuously to the tip. It then follows from 
(10.28) that dtdgfi extends continuously to the tip. A similar analysis shows 
that the same is true for dtdg/2- It then immediately follows from (10.28) 
that dfdgg a extends continuously to the tip. Differentiating (10.28) with 
respect to t one can repeat this analysis for all higher-order t-derivatives. 

An identical extension property is obviously true for P, Q, and all their 
derivatives. 

In order to prove (10.2), suppose that 

i-l 

dUaim = Y, a ^ t2j+1 +o(\tn, (10.43) 

3=0 

d k e9a {t^) = J2p a , j>k t 2 i + 0(\t\ 2i+1 ); (10.44) 

3=0 
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this clearly holds with i = 0. Then (10.28) gives 

i-1 

d t (td k e fa{t^)) = Y.aa,jM 2j+1 + 0(\t\ 2l+2 ), (10.45) 
3=0 

and by integration one recovers (10.43) with i replaced by i + 1. Inserting 
this into (10.28) one then finds 

i 

d t (d k e ga{tA)) = Y,^t 23+1 + 0(\t\ 2 ^ 2 ), (10.46) 

3=0 

and integration in t completes the induction step. The usual extension 
results finish the proof on C° Q (tp). Clearly the same results remain valid if v 
vanishes on an interval /, we simply note that all the estimates so far were 
uniform in 9 over /. □ 



11 Solutions with v < 1 and with uniform power- 
law 

Theorem 11.1 Let x be a solution of the Gowdy equations on Q(a,b,to) 
such that 

sup \tX t \ < 1 . (11.1) 

Q(a,b,t ) 

Suppose moreover that there exist positive constants e and Co such that on 
£l(a, b, to) we have 

\tX e \ < C \t\ e . (11.2) 

Then: 

(i) The solution belongs to the class %\ defined in Section 4, with a velocity 
function v = \v\\ and a position function (foo which are smooth except 
perhaps on the set 

d{9 : v{6) = 0} . 

(ii) There exists rj > such that for all initial data (x(to, •), X t (to, •)) 
satisfying 1 

\\(x(t ,-) - x(t ,-),X t (t ,-) ~ X t (t ,-))\\ H 3 @H 2 < K) 

the associated solution x also satisfies (11.1)-(11.2) (with perhaps dif- 
ferent values of e and C). Further the conclusions of Theorem 9.1 
apply. 
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Remark 11.2 The behavior of the solution near all points (0, ij)) with v(tp) = 
0, including {0} x (d{9 : v{9) = 0}), is described exhaustively in Theorem 10.1 
It is not clear whether or not the function ip^ of (3.24) is continuous at 
d{6 : v{6) = 0}. 

Remark 11.3 Every AVTD[ P ' 0) solut ion with v strictly smaller than one, 
and for which the error terms in (3.5)-(3.6) and in their derivative counter- 
parts (see (3.10)) are uniform in 9 satisfies (11.2), as well as (11.1) for some 
to close enough to 0. This follows immediately from (3.5)-(3.6). Further, all 
initial data satisfying the hypotheses of Theorem 8.6, or of [17, Theorems 1 
and 2], satisfy the hypotheses of Theorem 11.1. 

Proof: The arguments at the beginning of the proof of Theorem 8.6 show 
that v exists on [a, b] and is continuous there. We continue with a lemma: 

Lemma 11.4 Suppose that \tXg\ < C\t\ e for some e > and 

inf \tX t \ > 7 > , sup |£Xt| < 1 - 7 , (11.3) 

^KMl) n(o,6,tl) 

for some t±. Then Q(a,b,ti) can be covered by a finite number N of domains 
of dependence on which x is A VTD^'^ in appropriate coordinates on hyper- 
bolic space. (This implies smoothness ofv and of Voo; as those properties are 
invariant under isometries of the hyperbolic space.) The number N depends 
only upon e, Co and 7. 

We will actually prove a slightly stronger statement: 

Lemma 11.5 Under (11.3), suppose that there exists a sequence U such that 
the function 

2 2 

F(U) := sup Y.\ tk+lD o X ^°) + Y.\ tk+lD o X ^ e ) ( 1L4 ) 

ee[a+ti,b-ti] k=0 k=1 

goes to zero as ti goes to zero. Then the conclusion of Lemma 11. 4 holds, 
with the number N there depending upon 7 and the sequences {U}, {F(ti)}. 

Proof: We start by noting that under (11.2) we have F(ti) < C\U\ e by 
(10.8), so that Lemma 11.4 does indeed follow from Lemma 11.5 by setting 
U = 2-%. 

We wish to apply a result of Ringstrom [16, Theorem 9.1]. Set 

7=^min( inf \tX t \,l- sup \tX t \) . (11.5) 
8 fi(a,Mi) n(a,b,ti) 
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Let the constants e*,, k = 0, 1, 2 in [16, Theorem 9.1] be equal to 10, decreas- 
ing |ii| if necessary we can suppose that t± > —e~ T °, with To given by [16, 
Equation (9.1)]. Choose some tp £ (a, b), and define Q^VO to be a triangle 
with the top vertex at (0,tp) and slopes ±1/2: 

C^{ip) := {t>h, ip-2\t\ < 9 < ip + 2\t\} . (11.6) 

If ^ = a we set 

C£(V0 := {*>*!, Y> ~ \t\ < 6 < il> + 2\t\) , (11.7) 
while for ip = b we set 

CtM) -={t>ti,^- 2\t\ < e < v + \t\} . (ii.8) 

For < A < 1 let x^> A ) be defined as 

C£(V0 3 (t,0) -»• x^' A )(i,0) := x(At,V + A0) . (11.9) 

We have 

xf' A) (t,#) :=d t x^' x \t,0) =\X t (\t,iP + \6) , (11.10) 

similarly 

X^' X \t,6) := d e x^ x \t,6) = \X e (\t,iP + \6) , (11.11) 

with analogous formulae for higher derivatives. This, together with our 
hypothesis on the function F, shows that for A = Aj = U/ti we have for 

V> e (a, b) 

2 2 

sup ^i^x^Ktx.^+^I^^^Kti^) < n^i = ti)/l*il ■ 

ee[2ti ,-2ti] fe=0 fc=1 

(11.12) 

with similar results if ^ = a or ^ = 6. In particular at t = ti we have 

2 2 

sup ^ IDjX^Ktx,*) + \D k e X^ X) \(t a 9) -^ . (11.13) 

ee[-2\ tl \,2\ tl \] k=Q k=1 

Now, it is easy to show, using (5.6), that 

P 2 e + e 2P Q 2 w < \D e X t \ 2 + C\X t \ 2 \X e \ 2 , (11.14) 
Pee + e 2P Qi e < \D e X e \ 2 + C|X e | 4 . (11.15) 
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Similarly, 



Pm + e 2P Q 2 m < \D 2 e X t \ 2 + C(\X t \ 2 (\D e X e \ 2 + \X e \ 4 + \X e \ 2 ) 

+ \X e \ 2 \D e X t \ 2 ), (11.16) 
P 2 m + e 2P Q 2 eg < \D 2 X 9 \ 2 + C(\X \ 2 \D X t \ 2 + \X \ 6 + |A e | 4 ). (11.17) 

Equations (11.12)-(11.17) show that there exists Xj such that Equations (6.2) 
and (9.3) with k = 1,2, of [16] will be satisfied by the initial data for x(^> A ) 
for all A = A n < Xj with a multiplicative factor l/(2Ci) at the right-hand 
sides there, with a constant C\ to be made precise shortly. Decreasing Xj 
if necessary Ringstrom's energy e 2 corresponding to the current solution 
here will be smaller than l/(2Ci). Again decreasing Xj if necessary [16, 
Equation (9.4)] will hold with 27 there replaced by 47. It remains to satisfy 
Ringstrom's equation (9.3) with k = 0. Recall that the group of isometries 
of the hyperbolic space (J%, h) acts transitively on the unit tangent bundle 
of • This implies that there exists an isometry ip\ of {J%2 , h) such that 
the map ip\ox^^\ when written in the local coordinates (P, Q), will satisfy 

PM) = Q(t 2 , = d t QM) = , d t P(t 2 ,^) > 0. 

Now, Ringstrom's norms are not invariant under isometries. However, the 
objects appearing in (11.12) and (11.13) are, and those have been used to 
control Ringstrom's norms, so that the inequalities which have already been 
fulfilled still hold (in any case we could decrease A further to obtain the 
desired inequalities). It then follows by integration in 9, using the fact that 
the 9 derivatives of P and Q are already known to be small, that there 
exists d > such that Ringstrom's equation (9.3) with k = and with 
a multiplicative factor l/(2Ci) at the right-hand side there will hold on 
[tp — d + tj, tp + d — tj]. Decreasing d and Xj if necessary, we can extend the 
initial data from that last interval to smooth periodic initial data, without 
increasing all the relevant quantities by more than a factor C\. Applying 
Ringstrom's Theorem 9.1 the result follows for the map obtained by the 
evolution of the extended initial data. Uniqueness in domains of dependence 
establishes the claim on each — d,ip + d,tj). Reverting to the original 
ti, the result is obtained by covering Q(a,b,t±) by a finite number of sets 
Q(ip - d,ip + d,h). □ 



We note that an obvious modification of the argument just given estab- 
lishes the following version of Ringstrom's result [16, Theorem 9.1]: 
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Proposition 11.6 Let a < b, t < ; < 7 < 1, there exists 6(7) > such 
that if 

l<t P{to,-) < 1-7, 
2 2 
£ Itgi^Kto, •) + £ |*g^* t |(to, ■) < e , 

fc=0 fc=l 

on [a— \to\, b+\to\], then the solution is AVTD^'^ in£l(a,b,to), with velocity 
strictly positive and strictly smaller than one. □ 



We return to the proof of Theorem 11.1. The reader will note that the 
above proof has been worded to leave room for perturbing the data at t = t±, 
with ti as redefined in the paragraph following (11.5), while still satisfying 
Ringstrom's hypotheses; this is needed for the remainder of the argument. 

We set now ^ ^ 

7=-min(-,l- sup \tX t \) . 
8 6 n{a,b,t ) 

Since v is continuous, with tXt converging uniformly to v, the interval [a, b] 
can be covered by a finite number of intervals [a«, 6j] on which either 

7 < - inf \tX t \ , 

8 Q(ai,bi,ti) 

or on which ^ 

sup \tX t \ < — . 
n(ai,6i,ti) lz 

Let us call the latter intervals of type II, and the former of type I. Each 
of the [aj,6j]'s of type I can further be chosen to coincide with one of the 
intervals [ip — d, i/j + d] of the proof of Lemma 11.4 such that x is AVTD™ 
on 0,(ip — d,ip + d, £2)) where ti is the time given in the proof of Lemma 11.4. 
A sufficiently small pertubation of x at t = £2 leads again to an AVTD™ 
solution, thus satisfying (11.1)-(11.2) (with X there replaced by X, with 
possibly different constants C and e). The usual continuous dependence of 
solutions upon initial data on compact intervals of t shows that the same 
will remain true for sufficiently small perturbations of the initial data at 
t = t . 

On each interval of type II we have at t = t 2 

\tX e \ <Ct e < 1 -, 
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making ti smaller if necessary, then any sufficiently small perturbation of 
the Cauchy data at {£2} x [a% + £2, h — £2] leads to a solution x such that 

sup \tX t \ 2 + \tX e \ 2 < 2(1 + \) = \, 
n( ai ,bi,t2) 12 b 2 

where the factor 2 comes from Proposition 1.4. We can thus use Theorem 8.3 
to conclude x will satisfy (11.1)-(11.2). 

□ 



12 Existence and smoothness of v on an open dense 
set 

The proof of Theorem 1.3 will run in parallel with that of the following, 
more precise, statement: 

Theorem 12.1 Consider a solution x defined on Q(a,b,to) and set 

n := [sup \tX t \\ . 

There exists an open dense set Q, C [a, b] such that for every ip £ Q there 
exist a neighborhood ^Cfi of ip and an element G^ of the Geroch group, 
with the order of G^ less than or equal to n, such that G^x has on &^ a 
smooth velocity function < v = \v±\ < 1 and smooth position function Qoo- 

Proof: We start with a covering argument: We choose arbitrarily two 
points a < 9\ < 9 r < b, we set Iq = [61, r ] and we decompose £l(9i, 9 r , — \Io\/2) 
into a union of strips Qi, £l(9i, 9 r , —\Io\/2) = U^fii, where 

n t = {(t,x) g m,o r , -I/0I/2) I - ^ < t < -M} . 

Note that |/o|/2 is the height of Q(9i,9 r , —\Io\/2), and that the Q^s are 
pairwise disjoint. Let U = — and let Bi be the base line of Slj, 

Bi = {(t,x) efli\t = U, 9 t - \U\ <x<9 r + \U\}, i>l. 

We denote by Aio the left end point of Bi, and set 9\. = 9\ — |t,|, 9 n = 9 r + \ti\. 
We consider a partition £?i of Bi determined by the following sequence of 
points (see Figure 3) 

^ i = { ^iO , • • • , Aij , . . . , Ai^ mi } , 
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Figure 3: The points Ay and the triangles e^. 

where rrii = 2' + 2, Aij = (ti,9i t + ^rj)- Note that the length of each 

sub-interval of &i is 

Next, we decompose Qj into a union of triangles with 45 degrees slopes as 
follows: Let Cy denote the triangle with vertices at the points (Aij, Ai^j, Aij+i) 
defined above, where i > 1, j = 0, 1, . . . ,2 l + l. The set fij — UjJo 1 c «i ^ s then 
the union of 'upside down' triangles with vertices (Ai + i^j-2i A^, Ai + \^j), 
where i > 1, j = 1, . . . , 2* + 1. We denote those last triangles by dij. It 
follows that 

n, = (u?+^)|J(u?+%). 

Finally we relabel the Cy's and the dij's as e^: 

ej,2i+i := j = 0, . . . ,2* + 1 
e»,2j := J = 1, • • • ,2* + 1 

Lemma 12.2 Lei {fijj^Q 6e the decomposition offl(9i, 9 r , —\Iq\/2) described 
above. Let f be a nonnegative measurable function on fl(9i, 9 r , — | Iq |/2) with 

[ \t\f(t,9)d9dt < oo . 
Jn 

Then for any e > and jo £ N there exists j > jo such that fL contains a 
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set lj consisting of eight consecutive triangles eij with 

f f d9dt <e. 

Proof: Let J n \t\f(t, 9) dOdt = A. Using the decomposition £l(6i,0 r , -\I \/2) 



U^Qi we find 



oo „ 

A = ^2 \t\f{t,9)d9dt 

i=l 

>lIo \Y [ IMldedt. 



Thus, since the last sum is finite, there exists i(e) such that for all i > i(e) 
we have 

f f(t,6)d0dt<e 



for any given e. Let us show that there exist j and at least N = 8 consecutive 
triangles starting at e^ such that 

I f(t,0)d6dt<e. 

•' [J k=Q e i,i+k 

In fact we prove this assertion with any given N G N; we set 

en U e i2 ■ ■ ■ U eiN = c\ , 

e i,j U e «J+l • • • U e i,j+N-l = Cj , 
ei,m-N+l U ej im _Ar_|_2 . . . U ei m = C m -N+l , 

where m = 2 t+1 + 3. Let 



f f = Cj, [ f = Ej, 



then 



d + ... + C m . N+1 = E 1 + 2E 2 + ... + {N - l)£iv-i + N(E N + ... 
+ E m -N +3 ) + . . . + 2E m -i + E m 



<N(E 1 + ... + E m ) = N f f(t,6)d9dt < Nme. 
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Hence we obtain 

Ci + ■ ■ ■ + C m -N+i < Nme < ^ 
m-iV + 1 ~ m-N + 1 ~ 

which implies that there exists Cj such that Cj < Ne. □ 
We apply Lemma 12.2 to the function 

£ I 



f :=Y J \t k D k e Xe\ 2 + Y J \t k D k e X t \ 2 ■ 



k=0 k=l 

we are actually interested in I = 4, but the argument applies to any fixed 
£ € N, £ > 2. Point (iii) of Proposition 1.6 shows that / satisfies the 
hypotheses of Lemma 12.2, therefore there exists a sequence of domains 
uJi C fij such that 

/ > i~ >oo . 



/ 



The base b\ of UJi has length 4\U\, so that it can be written as 

bi = {U} x [^-2^1,^ + 2^1] , 

for some ipi £ [9i,6 r ]. We scale uj^ to a union of eight triangles with bottom 
edge lengths one, with the basis of the scaled set lying at t = —1, and (ipi, ti) 
mapped to (0, —1); we call Cj the resulting set, and we note that the top of 
Cj lies at t = — 1/2. We have 

/ / = / £ \t k DlX^ + £ |^X t ^| 2 , (12.1) 

with x^ uti \ etc., defined in (11.9). There exists an isometry of the hyper- 
bolic space into itself which maps x^'*^(0, —1) to the origin in the (P,Q) 
coordinate system. We apply this isometry to x^ l ' li \ and still use the same 
name for the resulting map. Since \t\ G [1/2, 1] on Co, isometry-invariance of 
(12.1) gives 

and the Sobolev inequality implies 

|| T t ^( 1 / , i;*i) || _j_ || T-i Y^itti) II n 

ll A 6> llc^- 2 (^) + IIM? A t llc*- 3 (u;) ^i->oo U . 
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Returning to the original Ui one thus has 

/i-2 1-2 



sup \t k+lD $Xe\ + E \ tk+lD e x t\ -i-oo • (12.2) 



fe=i / 
Let Ij = [tpi — 2\ti\,ipi + 2|tj|] and set 

Fe(Ii) ■= sup (j2 \t k+1 D k e X e \ + ^ |t* +1 I>JX t | > ) {ti,9) . 

ee[Vi-2|*i|,Vi+2|ti|] \fc=o fc=i / 

Equation (12.2) shows that i^O^i) approaches zero as ij tends to zero. We 
choose io large enough so that for all i > iq we have 

Clearly the same bound will then also hold for i^(7j) with < £ < 1. Since 
the sequence \tX t \(ti,ipi) is bounded, passing to a subsequence if necessary 
we can assume that there exists Voo such that \tXt\(ti,ipi) — > Vqo. Suppose, 
first, that < 1/200, then |iX t |(ij, ^j) < 1/100 for i large enough. By 
integration in 9 we have for 9 G ij 

||tX t | 2 (^,0)-|^t| 2 (^)| <±2 / i 2 |h(X tjJ D # X t )|cW<4sup|tX t |Fi(/ i ), 

(12.3) 

which goes to zero as i goes to infinity so that |iAi|(ij, 6) < 1/50 for 9 € ij. 
Now |t/i| < < 1/50, similarly for tf 2 , while |it/ 2 | < \tX e \ < F (Ii) < 

1/100, which proves that Theorem 8.3 applies, and shows that for all i > io 
(increasing io if necessary) the solution satisfies a power law decay in each 
of the — 2\ti\,tpi + 2\ti\,ti). The conclusions of Theorem 11.1 apply 
to show that the solution is AVTD^ > ' < ^ on each of the il(ipi — 2\ti\,ipi + 
2|tj|,tj), except perhaps a) for points in d{9 : v{9) = 0} at which Qoo 
may have discontinuities and/or v\ might fail to be smooth (even though 
it is continuous there) , or b) for a set of discontinuities of v\ introduced by 
applying back isometries to the isometry-transformed AVTD^'^- 1 solutions 
of Lemma 11.5. In any case Proposition 4.2 guarantees existence of an open 
dense subset of [9i,9 r ] with AVTD^'^ behavior there. 

Suppose, next, that 1/200 < < 1 — 1/200. A scaling argument as in 
the proof of Lemma 11.5 shows that, after applying a suitable isometry, the 
solution is AVTD^ >,< ^ in each of the £l(tfji — 2\ti\,ipi + 2|tj|,ij)'s, for i large 
enough, and hence again in an open subset of [9[, 9 r \. Applying the isometry 
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back to recover the original solution one obtains an open subset of [6i,6 r ] 
with AVTD^' Q) behavior. 

As the next possibility, consider the case in which 1 — 1/200 < Voo < 
1 + 1/200. Applying an isometry (pi we can assume that the map (pi o x, 
still denoted by x, satisfies Qt(U,ipi) = with P t {ti,ipi) - positive, so that 
1 - 1/100 < \t\P t (ti,ipi) < 1 + 1/100. Note that F e is invariant under 
isometries, so that F 2 (Ii) remains unchanged. Now, 

defa = d e {h{X t , e a )) = h(D e X t , e 2 ) + h(X t , D e e a ) , (12.4) 

and by integration, making use of (8.45) , one finds that 

-1/100 < \t\f 2 = \t\e p Q t (t h 8) < 1/100 

on Ii for i large enough, similarly 1 — 1/50 < \t\fi = \t\Pt(U,9) < 1 + 
1/50. Applying a Gowdy-to-Ernst transformation (3.20) (compare (8.33)) 
we obtain on Ii 

-1/50 < \t\P t = 1 - |t|P t < 1/50 , -1/100 < \t\ 9l = -\t\ 9l < 1/100 , 
as well as 

-1/100 < \t\g 2 = -\t\f 2 < 1/100 , -1/100 < \t\f 2 = -\t\g 2 < 1/100 . 

This shows that the hypotheses of Theorem 8.3 are satisfied by the initial 
data for i on I„ so that by Theorem 11.1 the map x will be AVTD^* Q) on 
an open dense subset of each of the Sl{ipi — 2\ti\,tpi + 2|tj|,tj)'s for i large 
enough. This gives an open subset of [0 h r ] with AVTD^' Q) behavior for 
x. To analyse the behavior of x we need first to perform back the Gowdy- 
to-Ernst transformation. Since 

P = -P-ln|t| , 

the existence and continuity properties of v are unchanged by the Gowdy- 
to-Ernst map. Next, we note that 

\Qt\ = \te p g 2 \ < CW 1 (12.5) 

for some e > 0, which by integration shows that the function Qoo of (3.6) 
exists and is a continuous function on each interval [ipi — \ti\,ipi + |tj|] with 
i large enough. Hence the map (P,Q) belongs to the class defined in 
Section 4. Now, the original map x is obtained from the one just described 
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by composing with an isometry of the hyperbolic plane, and is thus again 
in the class %\ by Proposition 4.2. One then obtains an AVTD^'^ map in 
a neighborhood of the set obtained by removing from [ipi — + |£j|] the 

countable set consisting of points where v has discontinuities, together with 
the boundary of the set where v = 1. 

We continue by induction: suppose that we have already established the 
claim for Voo < k + 1/200, and suppose that there exists k £ N such that 
k + 1/200 < < k + 1 — 1/200. Applying an isometry 0j to x, and still 
denoting by x the resulting map, we can assume that f 2 (U,ipi) = 0, and 
that A; + 1/50 < \t\P t < k + 1-1/50 on while \tX e \ + \tf 2 \ < 1/100 there. 
Applying a Gowdy-to-Ernst transformation we obtain 

—k + 1/50 < \t\P t < —{k — 1) — 1/50 on I h with \tf 2 \ < 1/100 . (12.6) 

It follows that \tX t \ < k on all ij's for i large enough. We note that (12.4) 
with a = 2 shows, by integration, that on I{ we have 

\t92\ = \tf 2 \ < CF^Ii) , while |tpi | = \t gi \ < F (Ii) holds trivially. 

(12.7) 

Next, by (5.7), 

DqXq = (d e gi - gl)ei + (d e g 2 + gih)^ 

= (-d e gi - /|)ei + (-d e f 2 + 9l f 2 )e 2 . (12.8) 

Similarly, 

D e X t = (defi - hhY\ + {deh + fm)e 2 

= (-d0fi-f 2 g 2 )e 1 + (-d e g 2 + f 1 f 2 -^)e 2 . (12.9) 

It is straightforward to check, using (12.7)-(12.9), that if one sets 

Fe(Ii) := sup (j2 \t k+1 D$X e \ + \t k+1 D$X t \) (U,0) , 

0efofc-2|ti|,Vi+2|ti|] \ fc=0 k=1 J 

then F 2 {Ii) — > as i — > oo. A similar calculation shows that the same is 
true for F 2 (Ii), and in fact for any higher order Fg if true for Fg. 

Returning to the proof of Theorem 12.1, it follows that the map x sat- 
isfies all the conditions needed to apply the result assumed to be true by 
the induction hypothesis. In particular the velocity function v and the po- 
sition function Qoo are smooth on each interval [ipi — \ti\,ipi + |tj|], with i 
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large enough, except perhaps for a countable number of points. Performing 
back the Gowdy-to-Ernst transformation one obtains a map with a velocity 
function v which has the same properties. Equation (12.5) becomes 



which shows, as before, existence of a continuous function on each in- 
terval [ipi — \U\,ipi + |tj|], with i large enough. Proposition 4.2 implies again 
the AVTDro'^ behavior in a neighborhood of on an appropriate subset of 
{0} x [ipi-\ti\,iPi + \ti\]. 

The possibility k + 1 - 1/200 < < k + 1 + 1/200 is handled similarly, 
and the induction step is completed. 

Since we have a uniform bound on \tX t \, the procedure stops in a finite 
number of steps. 

We have thus shown that every subinterval [9i,9 r ] of [a, b] contains an 
open set so that the solution is AVTD^'^ in a neighborhood thereof. The 
union of all such sets, as runs over all subintervals of [a, b], provides 

the desired open dense set Cl of Theorem 1.3. The set f2 of Theorem 12.1 can 
be taken to coincide with Cl. However, for the purpose of Proposition 12.3, 
one should keep in f2 those points at which the discontinuities in v arise 
from the action of the Geroch-group-inverse G^ 1 on Gx, when recovering 
x from Gx on — 2\ti\,ipi + 2|tj|,tj), by inverting the inductive method 
described above. □ 

The arguments given above together with Remark 1.2 can be used to 
obtain the following characterisation of points which are not in Q,: 

Proposition 12.3 Let £1 be the largest open set for which the conclusions 
of Theorem 12.1 hold. Then ip CI if and only if 

(i) either there exists e > such that for all t small enough 



(ii) or v exists and is continuous on an interval I containing ip, with ip 
belonging to the boundary of the set {v(9) = 0}. Further 



\Qt\ = \te p g 2 \ < Ce p < C\t\ k ^ 50 , 



(12.10) 



sup 

eeW-\t\^+\t\] 




2 




but Qoo is discontinuous at ip. 



□ 
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The key open question is that of existence of solutions for which the 
set of points exhibiting the properties described in Proposition 12.3 is not 
empty. The behavior described in (i) above seems to be a much more serious 
problem than the one in (ii). 

We are ready now to pass to the proof of the main result of our paper: 

PROOF of Theorem 1.1: Consider any point ip 6 [a, b\. The argument 
of the proof of Theorem 12.1 with ti = —2~ l , tpi = tp, shows the existence 
of a time < and of an element Gv, of the Geroch group such that 
G^x satisfies the hypotheses of Theorem 11.1 on fl(ip — \t^\,ip + \t^\,t^). 
(Note that the covering argument is not needed anymore in view of the 
hypothesis (1.6).) This proves (i) for Gx, and what has been said concerning 
the action of the Geroch group for the solutions under consideration proves 
(i) and (ii) for x. 

By compactness of [a, b] a finite covering by intervals (a{,bi) := (tpi — 
I^Vil'V' + \ttpi\) can be chosen, and (iii) readily follows. Point (iv) follows 
from the results in [4] (compare [6]). A small change of the initial data at 
t = to will lead to a small change of initial data at t = t^ i , hence to a small 
change of G^xtf^,-) and its derivatives at t^, and point (v) follows from 
Theorem 11.1. □ 
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